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ABSTRACT 


The  electron  density  distribution  through  ionized  clouds  produced 
by  Barium  releases  is  often  taken  to  be  Gaussian  for  want  of  direct 
evidence  to  the  contrary.  Multiple  frequency  hf  radar  observations 
(scattering  cross-section  and  Doppler  shift  versus  time)  have  been 
made  during  several  releases  for  the  primary  purpose  of  measuring 
peak  electron  density  as  a  function  of  time.  Observations  of  this 
sort  cannot  be  used  to  compute  the  electron  density  distribution 
directly  because  the  cloud  elongates  along  the  Earth's  magnetic  field 
and  thus  the  distribution  is  not  one  dimensional.  These  observations 
can  however  be  used  to  test  assumed  distributions  by  ray  tracing 
through  a  model  cloud  and  comparing  the  predicted  signatures  with  the 
observations.  A  detailed  comparison  of  this  nature  has  been  made 
using  a  time-dependent  Gaussian  cloud  and  the  ITSA  three-dimensional 
ray  tracing  program.  The  ray  tracing  results  have  been  compared  also 
to  that  of  a  hard  expanding  ellipsoid  model  to  study  the  effects  of 
the  ionization  lying  below  the  reflection  level.  Significant 
differences  have  been  noted  between  the  observed  and  the  synthesized 
signatures  from  the  ray  tracing.  The  differences  are  of  the  nature 
to  suggest  that  the  electron  density  variation  outward  from  the  center 
of  the  cloud  is  more  gradual  than  that  in  the  adopted  Gaussian  model. 
Defocussing  due  to  the  underlying  ionization  is  found  to  reduce  the 
scattering  cross-section  by  a  factor  that  increases  with  time  from 
3  at  the  peak  cross-section  to  30  just  prior  to  signal  drop-out. 
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SECTION  1.  INTRODUCTION 


The  generation  of  artificial  ion  clouds  in  the  upper  atmosphere 
has  proven  to  be  a  useful  technique  for  studying  plasma  interaction 
with  the  electric  and  magnetic  field  environment  of  the  Earth's  iono¬ 
sphere  and  magnetosphere.  The  technique,  originally  developed  at  the 
Max  Planck  Institute  in  Germany,  employs  rocket  releases  of  Barium 
vapor  which,  on  being  exposed  to  the  solar  ultraviolet  radiation, 
becomes  partially  ionized  (Haerendel  and  Lust,  1969) .  The  ionized 
and  the  non-ionized  components  of  the  Barium  vapor  form  into  two 
distinct  clouds,  an  ion  and  a  neutral  cloud.  A  number  of  radio  and 
optical  methods  have  been  devised  to  study  the  following  character¬ 
istics  of  the  clouds: 

a.  Size,  shape,  particle  number  densities  and  their  dis¬ 
tributions  as  a  function  of  time 

b.  Diffusion  and  drift  motions  including  the  dynamics  asso¬ 
ciated  with  the  gross  deformation  and  the  striated  struc¬ 
ture  exhibited  by  the  ion  cloud 

c.  Radio  wave  scattering  in  the  hf  and  vhf  regimes. 

The  knowledge  of  the  electron  (ion)  density  distribution  in  the 
ionized  cloud  is  of  fundamental  importance  for  a  proper  understanding 
of  the  ion  cloud  interaction  with  the  ionospheric  environment.  There 
have  not  been  as  yet  any  reliable  experimental  measurements  that 
would  directly  provide  this  information.  Some  evidence  has  been 
gathered,  however,  which  tends  to  suggest  a  Gaussian  profile  for  the 
distribution.  First,  the  Gaussian  model  has  stemmed  out  of  the  theo¬ 
retical  work  of  Holway  (1965)  which  is  based  on  the  ambipolar  diffu¬ 
sion  of  the  cloud  in  the  geomagnetic  field  neglecting  the  presence 
of  an  external  electric  field.  In  a  recent  study  by  Simon  (1970) ,  it 
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was  shown  that  taking  the  electric  field  into  account  would  still  lead 
to  a  Gaussian  distribution,  but  in  a  transformed  drift  frame  of  refer¬ 
ence.  The  experimental  evidence  that  seems  to  favor  the  Gaussian 
model  is  derived  from  the  optical  intensity  measurements  of.  the 
resonant  scattering  and  the  Thomson  scatter  radar  observations.  The 
deductions  from  optical  measurements  are  made  by  mierodens.itometer 
scans  of  the  filter  photographs  of  the  cloud.  The  variations  of  the 
line  integral  of  the  ions  along  selected  scans  are  found  to  fit  well 
to  Gaussian  distributions.  It  does  not,  however,  follow  immediately 
that  the  ion  density  distribution  itself  is  Gaussian.  Thomson  scatter 
radar  observations  provide  electron  density  profiles  along  the  antenna 
beam  direction  with  a  spatial  resolution  defined  by  the  pulse  length 
and  temporal  resolution  determined  by  the  integration  time.  The 
electron  density  profiles  measured  by  this  technique  using  the  AIO 
radar  during  one  of  the  SECEDE  I  tests  suggest  that  the  distribution 
is  at  least  approximately  Gaussian.  The  measurements  with  a  resolu¬ 
tion  of  2  km  are  rather  coarse,  however,  with  respect  to  the  cloud 
dimensions  and  are  subject  to  the  uncertainties  involved  due  to  pos¬ 
sible  energy  leakage  through  the  sidelobes.  Further  evidence  in  sup¬ 
port  of  the  Gaussian  distribution  has  come  from  an  analysis  by  Oetzel 
and  Chang  (1969)  who  found  that  the  observed  Faraday  fading  and  the 
history  of  the  center  point  electron  density  derived  from  hf  measure¬ 
ments  are  consistent  to  a  greater  extent  with  the  Gaussian  model  than 
with  other  distributions. 

In  spite  of  all  the  positive  evidence  in  support  of  the  Gaussian 
model,  the  true  distribution  is  still  a  matter  of  speculation  because 
of  an  unexpected  behavior  of  the  hf  Doppler.  It  is  expected  that  an 
expanding  Gaussian  cloud  will  lead  to  a  Doppler  trace  that  goes  neg¬ 
ative  shortly  before  the  signal  dropout  on  a  given  operating  frequency. 
The  observed  Doppler  records,  however,  exhibit  no  such  negative  tail. 


For  this  to  occur,  two  possibilities  suggest  themselves  immediately: 
one,  of  course,  is  that  the  distribution  is  unlike  Gaussian,  and  the 
other,  which  seems  more  probable,  is  that  the  signal  is  too  weak  to 
be  detected  during  the  contracting  phase  of  a  given  plasma  frequency 
contour,  thus  washing  off  the  expected  negative  tail. 

The  problem  taken  up  in  this  report  is  to  investigate  the  two 
possibilities  by  synthesizing  the  radar  signatures  both  in  Doppler  and 
signal  strength  by  means  of  a  three  dimensional  ray  tracing  program 
and  Gaussian  modeled  ion  clouds.  The  results  will  show  to  what  ex¬ 
tent  the  sensitivity  of  the  receiving  system  limits  the  Doppler  trace, 
and  how  accurately  a  Gaussian  model  could  simulate  an  ion  cloud.  The 
3-D  ray  tracing  program  used  here  was  developed  at  ESSA  by  Dr.  R.M. 
Jones  and  has  proven  to  be  quite  versatile  in  its  scope.  A  summary 
description  of  this  program  is  given  in  Section  2.  Section  3  deals 
with  the  ion  cloud  modeling  in  the  geomagnetic  dipole  coordinate  sys¬ 
tem  as  desired  by  the  ray  tracing  program.  The  results  of  the  ray 
tracing  are  presented  in  Section  4  along  with  those  derived  for  the 
case  of  an  expanding  hard  ellipsoid  model.  Experimental  results  are 
described  in  Section  5»  Section  6  is  devoted  to  a  comparison  between 
the  theory  and  observations  with  the  summary  and  conclusions  given  in 
Section  7. 


SECTION  2.  THREE  DIMENSIONAL  RAY  TRACING  PROGRAM 


2.1  INTRODUCTION 

A  versatile  computer  program  for  tracing  rays  through  a  medium 
whose  index  of  refraction  varies  in  three  dimensions  has  been  de¬ 
veloped  by  R.  M.  Jones  at  ESSA  Laboratories  in  Boulder.  The  program 
was  intended  for  ray  tracing  applications  in  the  ionosphere  and  could 
successfully  simulate  a  variety  of  radio  propagation  experiments.  It 
was  written  in  FORTRAN  language  and  was.  made  available  on  request  to 
prospective  users  in  two  versions;  one  to  run  on  a  CDC  3600  and  the 
other  on  an  IBM  7090.  A  detailed  description  of  the  program  along 
with  instructions  that  include  annotated  listings  of  all  routines, 
deck  set-ups  and  the  form  of  input  and  output  has  been  published  in 
the  ESSA  reports  by  Jones  (1966,  1968) .  The  efficiency  of  the  pro¬ 
gram  in  computing  the  rays  has  been  evaluated  by  Lemanski  (1968) . 

2.2  GOVERNING  EQUATIONS 

The  equations  that  govern  ray  paths  in  a  three  dimensional  ani¬ 
sotropic  medium  have  been  derived  using  the  method  of  Hamilton  by 
Haselgrove  (1954) .  The  equations  were  expressed  in  a  form  suitable 
for  integration  by  standard  numerical  methods  using  a  digital  com¬ 
puter.  The  Jones  program  calculates  the  ray  paths  by  integrating 
the  following  six  differential  equations  which  are  a  slight  modifi¬ 
cation  to  those  described  by  Haselgrove. 
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'The  variables  r,  0,  and  0  are  the  spherical  polar  coordinates  of  a 

point  on  the  ray  path;  V  ,  V ,  and  V  are  the  components  of  the  wave 
c  r  o  0 

normal  in  the  r,  6,  and  0  directions,  normalized  so  that 


V2  +  V2.  +  V2  =  Real  (n2)  (2-7) 

r  6  0 

where  n  and  n'  are  the  complex  phase  and  group  refractive  indices 
which  are  based  on  Apple ton -Har tree  formula.  The  independent  vari¬ 
able  t  is  group  path  length  in  these  equations  whereas  it  is  phase 
path  length  in  Haselgrove  equations.  This  change  of  independent 
variable  has  the  advantage  of  speeding  up  the  program  by  automat¬ 
ically  decreasing  the  step  length  in  real  path  near  the  reflection 
height.  In  addition  to  the  six  basic  equations  necessary  to  cal¬ 
culate  the  ray  path,  the  program  integrates  the  following  to  give 
the  phase  path,  the  absorption,  and  the  Doppler  shift,  respectively. 
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2 .3  DESCRIPTION  OF  THE  PROGRAM 

The  program  traces  the  path  of  a  radio  wave  through  the  ionosphere 
when  given  along  with  the  ionospheric  models,  the  transmitter  coor¬ 
dinates,  the  operating  frequency,  the  direction  of  transmission,  the 
height  of  the  receiver  and  the  maximum  number  of  hops  desired.  The 
general  set  up  of  the  program  can  take  both  the  electron  collisions 
and  the  Earth' s  magnetic  field  into  account  and  trace  the  two  types 
of  ray  paths  (ordinary  and  extraordinary) ;  It  can,  at  our  option, 
drop  either  or  both  of  the  collisions  and  the  magnetic  field.  The 
automatic  homing  feature,  common  to  some  ray  tracing  programs,  is 
not  built  into  this  program  and  it  had  to  be  achieved  by  trial  and 
error,  knowing  the  starting  point  for  the  approximate  direction  of 
transmission  to  the  receiver. 

The  distributions  of  electron  density,  collision  frequency  and 
Earth's  magnetic  field  constitute  the  input  models  that  describe  the 
ionosphere.  Routines  for  several  ionospheric  models  have  been  made 
available  with  the  program.  Five  analytic  models  and  tabular  or 
true-height  analysis  profiles  constitute  the  supplied  electron  den¬ 
sity  routines.  Two  collision  frequency  routines  exist;  one  for 
tabular  profiles  and  the  other  for  a  constant  collision  frequency. 
There  are  two  routines  that  model  the  Earth's  magnetic  field;  one 
is  for  an  Earth-centered  dipole,  and  the  other  is  based  on  constant 
dip  and  gyro -frequency.  The  model  routines  are  arranged  in  the 
programs  in  such  a  fashion  that  they  can  be  replaced  with  any  de¬ 
sired  models  with  utmost  ease  provided,  of  course,  the  new  routines 
follow  the  format  already  set. 


The  numerical  integration  of  the  ray  tracing  equations  can  be 
carried  out  by  using  either  one  of  the  two  routines  available;  one 
based  on  Runge-Kutta  method,  and  the  other  on  Adams-Moulton  method. 


The  integration  routine  by  Adams -Moulton  method  has  a  built-in  mech¬ 
anism  for  checking  errors  and  adjusting  the  integration  step  size 
accordingly.  The  maximum  allowable  error  in  any  single  step  for  any 
of  the  equations  integrated  in  the  routine  is  specified  in  the  form 
of  an  input  parameter.  The  step  size  is  decreased  to  gain  required 
accuracy  if  the  error  is  larger  than  the  specified  maximum  and  in 

creased  to  reduce  the  computing  time  if  the  error  is  smaller.  A 
-4 

value  of  10  is  commonly  used  for  the  maximum  allowable  error  as  a 
compromise  between  cost  and  accuracy. 

The  program  uses  two  spherical  polar  coordinate  systems;  one  is 
geographic,  and  the  other  computational.  The  input  data  such  as  the 
location  of  the  transmitter  and  the  geomagnetic  north  pole  is  given 
in  geographic  coordinates;  the  transformation  to  the  computational 
system  is  carried  out  within  the  program.  When  the  dipole  model  is 
used  for  the  Earth's  magnetic  field,  the  computational  system  is  a 
geomagnetic  coordinate  system  and  both  electron  density  and  colli¬ 
sion  frequency  must  be  defined  in  geomagnetic  coordinates. 

The  input  data  to  the  program  is  arranged  in  a  one  dimensional 
array  (W  array)  and  read  in  with  one  parameter  on  each  card.  This 
array  provides  all  the  information  needed  for  the  ray  trace  includ¬ 
ing  the  parameters  associated  with  the  analytic  models  for  electron 
density,  collision  frequency  and  magnetic  field.  The  output  is 
available  in  various  forms  individually  or  in  any  combination. 

During  the  course  of  the  ray  trace  the  program  will  print  information 
giving  the  position  of  the  current  ray  path  point,  the  wave  normal 
direction  and  cummulative  values  of  the  quantities  like  group  path, 
phase  path,  absorption,  and  Doppler  shift.  The  main  results  of  the 
program  can  be  obtained,  if  desired,  in  the  form  of  punched  cards 
called  ray-sets  which  can  be  used  as  input  for  other  programs.  The 
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ray  paths  can  be  plotted  in  two  projection  planes?  a  ground  plane 
and  any  given  vertical  plane.  It  is  possible  to  show  the  lateral 
deviations  in  a  blown-up  scale  on  the  ground  plane  plot. 

2.4  APPLICATION.  TO  PLASMA  CLOUDS 

The  ray  tracing  program  has  been  used  here  to  investigate  the 
scattering  behavior  of  an  expanding  Gaussian  plasma  cloud  by  synthe¬ 
sizing  hf  direct  echo  signatures  both  in  Doppler  shift  and  signal 
strength.  The  rays  that  have  been  traced  are  intended  to  simulate 
the  signatures  of  test  Mulberry  (Pre-SECEDE)  for  which  high  quality 
Doppler  data  is  available.  The  ultimate  objective  here  is  to  compare 
the  results  with  observation  and  thus  test  the  validity  of  the 
Gaussian  model  for  the  ion  cloud.  The  program  has  been  used  in  the 
mode  that  takes  the  collisions  and  magnetic  field  into  account  as 
well  as  in  the  modes  that  neglect  either  one  or  both.  The  ray  paths 
have  been  calculated  for  all  the  three  types  of  rays,  namely,  ordin¬ 
ary,  extraordinary  and  no-field. 

The  models  adopted  for  the  electron  density  of  the  cloud,  the 
collision  frequency  and  the  magnetic  field  are  all  analytic  through¬ 
out.  The  ion  cloud  model,  discussed  in  detail  in  the  next  section, 
is  represented  by  a  Gaussian  distribution  cylindrically  symmetric 
about  the  direction  of  Earth's  magnetic  field  and  governed  by  aniso¬ 
tropic  diffusive  expansion.  The  relevant  Gaussian  parameters  that 
determine  the  electron  density  in  time  and  space  are  based  upon  ex¬ 
perimental  measurements.  The  electron  collisions  are  represented  by 
a  piece-wise  exponential  function  with  parameters  that  best  fit  the 
models  published  by  Thrane  and  Piggot  (1966).  It  turned  out,  however, 
that  the  collisional  effects  are  negligible  and  never  amount  to  more 
than  a  fraction  of  a  dB  in  absorption.  This  is  due  to  the  fact  that 
the  ion  cloud  considered  here  was  released  in  the  F  region  where 
collisions  are  small  and  at  twilight  when  the  contribution  from  the 


ambient  ionosphere  is  still  negligible.  The  dipole  model  is  chosen 
for  the  Earth's  magnetic  field  with  the  geomagnetic  north  pole  lo¬ 
cated  at  78.3°N  and  291°E  (geographic).  The  equatorial  value  of  the 
gyrofrequency  is  taken  to  be  0.8  MHz  and  this  together  with  the  di¬ 
pole  model  determines  its  value  at  any  point  in  space. 

The  method  of  Adams-Moulton  with  relative  error  check  was  chosen 

out  of  the  four  options  provided  in  the  program  for  integrating  the 

ray  trace  equations.  The  maximum  allowable  single  step  error  was 

-4  -6 

given  a  value  in  the  range  10  to  10  depending  on  how  stringent 

the  ray  tracing  conditions  are  at  a  given  time  of  the  cloud's  life. 

—6  -6 

It  was  assigned  a  smaller  value  (5  x  10  or  10  },  to  keep  uniform 

accuracy  as  judged  frorr  the  first  column  of  the  printout  (see  Appen¬ 
dix  A) ,  during  the  late  phase  when  the  plasma  frequency  contour  sup¬ 
porting  the  signal  return  is  small  in  size  and  rapidly  collapsing 
to  the  center.  Whenever  this  parameter  was  altered,  suitable  changes 
have  also  been  made  in  other  related  parameters  like  initial  integra¬ 
tion  step  size  and  its  minimum  limit. 

The  ray  tracing  results  were  printed  out  in  all  cases,  but  plots 
in  the  two  possible  planes,  vertical  and  ground,  have  been  made  only 
for  some  selected  cases  of  interest.  A  sample  plot  and  a  printout 
are  presented  in  Appendix  A.  The  supplementary  information  on 
polarization,  phase  path  and  group  path,  obtained  along  with  the 
Doppler  results,  will  not  be  discussed  here  as  they  fall  beyond  the 
scope  of  this  report. 
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SECTION  3.  PLASMA  CLOUD  MODEL 


3.1  AMBIPOLAR  DIFFUSION  MODEL 

The  ion  clouds  of  interest  here  are  generated  in  the  F  region 
by  means  of  rocket  releases  of  metallic  Earium  vapor.  The  cloud,  on 
formation,  is  subject  as  a  whole  to  electrodynamic  drift  due  to  the 
presence  of  external  electric  field  while  at  the  same  time  expanding 
in  a  manner  determined  by  diffusion.  The  charge  neutrality,  that  is 
n^  =  n^,  is  satisfied  in  the  cloud  since  the  Debye  shielding  length 
is  much  smaller  than  the  typical  cloud  dimensions  and  the  diffusion 
is  of  ambipolar  nature.  The  problem  of  ambipolar  diffusion  of  the 
cloud  in  the  geomagnetic  field,  neglecting  the  external  electric  field, 
has  been  studied  by  Holway  (1965) .  The  diffusion  equation  has  been 
solved  in  cylindrical  coordinate  system  (r,  0,  z)  with  z  axis  taken 
along  the  Earth's  magnetic  field.  In  this  system,  because  of  the 
cylindrical  symmetry,  the  azimuthal  component  of  the  particle  flux  is 
non-divergent  and  it  is,  therefore,  necessary  to  consider  only  the 
components  in  the  z  and  r  directions.  Using  the  above  criterian  and 
a  relation  that  the  ion  flux  is  the  same  as  the  electron  flux  in  both 
the  radial  and  logitudinal  directions,  Holway  derives  the  equation  for 
the  ambipolar  diffusion  of  the  ionized  cloud  as: 


dn  _  D||  1  d  dn 
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where  n  is  the  ion  (electron)  number 
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The  solution  to  Equation  3-1  for  an  ionized  release  of  ions,  which 
is  initially  spherical  with  a  Gaussian  distribution  of  initial  effec¬ 
tive  radius  h  ,  is: 

o 


n  (r,  z,  t)  = 


-yT~— 2I  axp  [-  <z/zi)2  ‘  (r/ri>2]  (3-2) 
L*  Vi  J 


where 


z  =  (h2  +  4D. .  t)1^2  and  r,  =  (h2  +  4D  .  t/Q.  are  the 

1  o  11  1  o  1 1  le 


longitudinal  and  transverse  Gaussian  radii. 

An  extension  has  been  made  recently  to  the  above  solution  by 
Simon  (1970)  by  taking  the  presence  of  an  external  electric  field  in¬ 
to  consideration.  The  diffusion  equation  in  this  case  has  been  form¬ 
ulated,  after  making  considerable  simplification  by  Simon  as: 

|f  =  V±  V2±  N  +  P(|  y2  N  +  a  Uq  .VN  +  0  UQ  x  b  .VN  (3-3) 

The  solution  to  Equation  3-3,  assuming  that  the  plasma  is  created  at 

the  origin  at  t  =  0  and  the  total  number*  of  ions  created  at  that  time 

is  S  ,  is 
o 


N  (x,  t)  = 
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Note  that  Equations  3-3  and  3-4  retain  Simon's  notation.  If  the 
initial  condition  is  changed  from  a  5  function  to  that  of  a  spherical 
Gaussian  distribution  with  radius  h^,  the  Gaussian  radii  of  Equation 
3-4  would  transform  to  the  same  form  as  of  Equation  3-2  in  the  manner 
discussed  by  Holway  (1965).  The  solution  given  in  Equation  3-4  is  a 
cylindrically  symmetric  Gaussian  moving  with  velocity  (-  -  0u  x  b)  . 

Since  the  Doppler  associated  with  uniform  drift  velocity  can  be  taken 
into  account  explicitly  outside  the  ray  tracing  program.  Equation  3-2 
has  been  adopted  here  for  the  electron  density  distribution.  The 
origin  about  which  the  expansion  takes  place  was  allowed  to  shift, 
however,  on  the  basis  of  the  trajectory  information  obtained  from  the 
optical  measurements. 

3.2  MODEL  TRANSFORMATION  TO  GEOMAGNETIC  COORDINATE  SYSTEM 

The  ray  tracing  program  calls  for  electron  densities  and  their 

gradients  in  terms  of  geomagnetic  spherical  dipole  coordinates.  This 

requires  a  correspondence  between  the  cloud  centered  cylindrical 

coordinate  system  and  the  Earth -centered  dipole  coordinate  system. 

The  transformation  can  be  made  first  by  relating  two  Cartesian  systems, 

one  cloud  centered  (x,  y,  z)  and  the  other  Earth  centered 

(x^,  y^,  z^) ,  since  they  are  directly  related  to  respective  polar 

systems.  The  system  (x,  y,  z)  is  defined  such  that  z  is  anti -parallel 

to  the  magnetic  field  direction  and  y  is  directed  toward  geomagnetic 

east  in  the  northern  hemisphere.  In  (x  ,  y  ,  z  ) ,  z  is  along  the 

m  m  m  m 

dipole  axis  and  directed  toward  the  north  pole,  (x  ,z  )  defines  a 
o  mm 

reference  (0  )  geomagnetic  meridian  plane  and  y^  is  toward  geomagnetic 

east.  In  general,  it  involves  one  translation  and  two  rotations  to 

bring  correspondence  between  the  two  systems,  as  shown  in  Figure  3-1. 

The  translation  is  along  the  vector  joining  the  two  origins  and  the 

rotations  are  in  two  orthogonal,  geomagnetic  latitude  and  meridian, 

planes.  Let  the  cloud  center  in  geomagnetic  coordinates  be 
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magnetic 

FIELD 

LINE 


Figure  3-1  Illustration  to  Show  the  Coordinate  Transformation  from 
Earth  Centered  Geomagnetic  Dipole  System  (x  ,y  ,z  )  to 
Cloud  Centered  System  (x,y,z).  The  Correspondence  Be¬ 
tween  the  Two  Systems  Requires  One  Translation  and  Two 
Rotations. 


(R  ,  0  ,6  ),  9  being  geomagnetic  co-latitude.  When  the  system 

o  mo  mo  mo  ~ 

(x  ,  y  ,  z  )  is  translated  to  the  cloud  center  and  rotated  by  rf> 
m  m  m  J  vmo 

about  the  Z  axis,  an  intermediate  system  (x' ,  y>  ,  z>)  with  y'  co¬ 
rn  m  m  m  m 

inciding  with  y,  would  result  as  shown  in  Figure  3-1.  (x' ,  y' ,  z') 

m  m  m 

is  related  to  (x  ,  y  ,  z  )  as: 

m  m  m 


cos  <{>  sin  4>  0 

mo  mo 


X“1 


■sin  <j>  cos  4>  0  v 

mo  mo  n 


(3-5) 


The  matrix  elements  a',  b' ,  and  c'  are  the  result  of  the  translation. 

Now  a  rotation  by  an  angle  v  ,  which  is  equal  to  (90-1+0  )  as  shown 

J  ^  Amo  ^  mo 

in  Figure  3-1,  in  the  plane  (x' ,  z' )  about  the  y >  axis  would  bring  the 

mm  m 

intermediate  system  to  coincide  with  (x,  y,  z) .  The  relation  between 
the  two  systems  is  given  as: 


cos  x. 


sm  v 


0  -sin  x 


0  cos  X 


(3-6) 


combining  Equations  3-5  and  3-6  we  have 


cos  4>  cos 
mo 


sin 


sin  4>  cos 
mo 


cos  4> 


X  -sin  v  x 
^no  Amo  m 


cos  4>  sin  v  sin  4>  sin  v  cos  v  z 
mo  Amo  mo  Amo  AmoJ  l 


(3-7) 
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The  elements  a,  h  and  c  denote  the  translation  with  respect  to  x,  y, 
and  z  and  can  be  evaluated  by  considering  the  origin  of  the  cloud 
where: 


x 

y  = 


X  = 

x  = 

(R 

+ 

h) 

sin 

0 

cos 

m 

mo 

© 

mo 

y  = 

y  = 

(R_ 

h) 

sin 

0 

sin 

m 

mo 

© 

mo 

z  = 

z  = 

(R_ 

+ 

h) 

cos 

e 

m 

mo 

© 

mo 

:  the 

Earth, 

h  : 

Ls 

the 

i  height 

and 

(3-8) 


©  mo  rmo 

geomagnetic  colatitude  and  longitude  of  the  cloud  center.  Substituting 
Equation  3-8  into  Equation  3-7  gives 


a  =  -x  cos  <f>  cos  v  -  y  sin  <{>  cos  y  +  2  sin  X 

mo  mo  Amo  mo  mo  mo  mo  mo 


b  =  x  sin  4>  -  y  cos  <f> 

mo  mo  mo  mo 


(3-9a) 
(3  -9b) 


c  =  x  cos  sin  y  -  y  sin  4>  sin  y  -  z  cos  > 

mo  mo  Amo  mo  mo  "Tno  mo  mo 

(3 -9c) 


The  position  of  the  cloud  center  is  made  available  usually  in  geo¬ 
graphic  coordinates  and  they  are  converted  to  geomagnetic  coordinates 
using  the  following  expressions: 


cos  0  =  sin  0  sin  0  +  cos  0  cos  0  cos  (<j>  -<j>  ) 

mo  on  on  o  n 


(3-10a) 


sin  $ 


cos  0  sin  (<}>  -  <j>  ) 

_ 0 _ Q  n 

sin  0 


(3-10b) 
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Where  (0  ,  d>  )  are  the  geographic  latitude  and  longitude  of  the  cloud 
o  o  ~ 

center  and  (6  ,  $  )  are  the  geographic  latitude  and  longitude  of  the 
n  n 

geomagnetic  north  pole  (note  that  0^q  is  the  geomagnetic  colatitude) . 


When  the  ray  tracing  program  refers  to  the  electron  density  sub¬ 
routine  with  geomagnetic  coordinates  (R, 9  ),  they  are  transformed 

m  m 

to  corresponding  (x,  y,  z)  through  (x^,  y  ,  z^) .  Further,  by  con¬ 
verting  the  cylindrical  coordinates  of  Equation  3-2  to  corresponding 
(x,  y,  z)  we  have: 


n(R,0  ,4>  )  =  n  (x,y,z)  = 
mm 


N. 

L  7r3/2 


exp. 


r  f  2  .  2  v  2  -i 

-lx  tv)-.  2_ . _ 


(3-11) 


3.2.1  DENSITY  GRADIENT  FIELD 

The  electron  density  partial  derivatives  are  evaluated  by  using 
the  following  matrix  relations. 


9n 

9x 

m 

9y 

m . 

9z 

m 

9n 

9  R 

9  R 

9R 

9R 

9x 

m 

9n 

9x 

m 

9y 

Jm 

9z 

m 

9n 

90 

m 

90 

m 

90 

m 

90 

m 

9y 

m 

9n_ 

9x 
.  m 

aY 

m  . 

9z 

m 

9n 

94> 

m  _ 

9cf> 

L  m 

9<j> 

m 

9<j> 

m 

9z 

L  m 
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sin  9  cos  0 


sin  9  sin  0  cos  9 

m  rm  n 


R  cos  9  cos  0  R  cos  9  sin  <4  -R  sin  & 
m  m  m  n 


-R  sin  9  sin  <4 
m  ^m 


R  sin  0  cos  * 
m  ^m 


(3-13) 


'  d% _ 

dy 

dz 

3x 

dx 

dx 

m 

m 

m 

9x 

hL- 

dz 

3y 

Jm 

*Ym 

m 

aym 

m 

9x 

£y_ 

9z 

dz 

9z 

9z 

-g-Y... 

dz 

m 


(3-14) 


cos  0mQ  cos  Xmo  -sin  0mQ  cos  0  sin  v 


s-*-n  ^tnn  cos  Xmri  cos  0  sin  0  sin  v 
mo  mo  rmn  'mn 


l-sm  v 


cos  x 


(3-15) 


n (R, 9  ,<j>  )  [~-  2x/r  2  1 
m  m  L  1  J 


n(R,0  ,4>  )  f-  2y/r  2  "1 
m  m  L  1  J 


(3-16) 


n (R/ 0  t$  )  [*-  2z/z  2  1 
m  m  L  1  J 


The  time  derivative  3n/3t  is  evaluated  by  writing  n  as 


N.  (i-e"t/T 
n(x,y,z,t)  =  ~^2  "S  2 


2  2 
x  +  v 


^..v—  2 

2  r 

z  r  1 

1  1 


e  z. 


f  (t) 


g(t) 


H  ( t) 


(3-17) 


j, _  ran. 

3/2  |_3f 


3  f  +  3_n_  3g_  ,  3n  3H 
3t  3g  3t  3H  3t 


(3-18) 


“at  +  fH  ft-  +  fg  ft" 


-t/T 


r  z  r 
1  1 


-TT  <4  zi  D1  +  2DII  rl2> 
Z1  ri 


(3 -19a) 


i3L_  _ 
3 1 


4Di9  [(x2  +  y2)/ri4] 


(3 -19b) 


4  D j  |  H  (z  /z14) 


(3 -19c) 
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3.3  ION  CLOUD  INPUT  PARAMETERS 


The  coordinates  of  the  center  of  the  ion  cloud  and  the  parameters, 
total  ion  content  N^,  ionization  time  constant  t,  initial  Gaussian 
radius  hQ,  and  longitudinal  and  transverse  diffusion  coefficients 
Dj j  and  form  the  input  to  set  up  the  electron  density  distribution 
for  the  ray  tracing.  The  numerical  values  of  the  parameters  are 
appropriate  to  test  Mulberry  of  Pre -SECEDE  and  are  taken  from  Minkoff 
(1970) .  The  positional  data  obtained  by  photographic  means  was 
available  in  the  form  of  look  angles  (azimuth  and  elevation)  from 
different  observation  sites.  These  observations  taken  from  pairs  of 
stations  have  been  used  to  determine  the  coordinates  of  the  cloud 
following  a  vector  method  described  in  Appendix  B.  The  trajectory 
information  is  shown  in  Figure  3-2  for  the  period  of  interest  of  ray 
tracing. 

The  total  ion  content  N^  is  obtained  by  integrating  Equation  3-2 
over  all  space,  from  which: 

N.  =  n  (t)  tt3/2  (h2  +  4D  t)1/2  (h2  +  4D.  t)  (3-20) 

i  p  o  ||  o  1 

The  peak  electron  (ion)  density  n  (t)  is  taken  from  hf  radar  measure- 

P 

2  1/2  2  1/2 

ments.  The  Gaussian  radii  (h  +  4D  t)  and  (h  +  4D,  t)  are 

o  U  o  _L 

determined  by  means  of  microdensitometer  scans  of  the  filter  photo¬ 
graphs  of  the  cloud.  Figure  3-2  also  presents  the  peak  electron 
density  variation  with  time  and  the  results  from  the  densitometer- 
scans.  The  numerical  values  obtained  for  the  parameters  are  given 
in  Table  3-1. 
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ALTITUOEIKM 


Figure  3-2  The  Trajectory,  the  Peak  Electron  Density  and  the  Results 
From  Microdensitometer  Scans  of  the  Filter  Photographs  of 
the  Cloud.  (Filled  and  Open  Circles  in  the  Trajectory  Data 
Represent  Determinations  from  Two  Pairs  of  Optical  Stations) 


TABLE  3-1 


NUMERICAL  VALUES  FOR  ION  CLOUD  PARAMETERS 


Test 

h  (KM) 
o 

D  (KM2/SEC) 

2 

(KM  /SEC) 

N. 

l 

Mulberry 

3.0 

0.13 

0.045 

3.0  x  1024 

The  reader  is  referred  to  Minkoff  (1970)  for  details  on  how  these 
parameters  are  deduced. 
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SECTION  4.  SYNTHESIS  OF  RADAR  SIGNATURES 


Two  methods  have  been  used  for  synthesizing  radar  signatures  for 
the  test  Mulberry;  one  involves  the  rigorous  ray  tracing  technique, 
and  the  other,  an  approximate  technique  meant  primarily  for  the  purpose 
of  comparison,  is  based  on  an  expanding  hard  ellipsoid  model.  This 
section  is  devoted  to  a  description  of  the  two  methods  and  to  a  dis¬ 
cussion  of  the  results  obtained. 

4.1  RAY  TRACING  METHOD 

There  are  two  steps  involved  in  synthesizing  a  radar  signature. 

The  first  step  is  to  determine  the  ray  that  homes  in  on  the  receiver 
when  the  transmitter-cloud-receiver  configuration  is  given.  The 
Doppler  shift  of  interest  is  the  one  associated  with  this  ray  and  it 
is  listed  in  the  output  of  the  ray  tracing  program  as  a  cumulative 
sum  along  the  ray  path.  The  next  step  is  to  determine  the  signal 
strength  associated  with  this  ray.  It  is  evaluated  in  terms  of  radar 
scattering  cross  section  of  the  cloud  since  it  is  a  convenient  pa¬ 
rameter  to  compare  different  models  and  to  gain  some  insight  into  the 
scattering  behavior  of  the  ion  clouds.  The  process  of  homing- in  on 
the  receiver  proceeds  in  steps  and  is  as  follows:  First,  the  eleva¬ 
tion  and  azimuth  of  the  cloud  center  is  calculated  with  respect  to 
the  transmitter  location.  This  provides  the  initial  guess  for  the 
direction  of  transmission  at  the  starting  time.  Using  this  as 
reference  direction,  a  matrix  of  nine  rays  (3x3  in  elevation  and 
azimuth)  defining  a  rectangular  solid  angle  are  shot  normally  as  a 
first  trial.  The  angular  increments  in  elevation  and  azimuth  are 
selected  such  that  the  impact  points  form  a  contour  that  encloses  the 
receiver  location.  The  points,  in  general,  describe  a  skewed  quad¬ 
rilateral  because  of  the  nonlinear  nature  of  the  cloud  defocussing. 
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The  next  trial  uses,  in  general,  four  rays  for  which  the  elevation  and 
azimuth  angles  are  arrived  at  by  interpolating  between  the  constant 
elevation  and  constant  azimuth  lines  associated  with  the  impact  points. 
The  rays,  in  this  trial,  converge  toward  the  receiver  to  come  usually 
within  25  km.  Proceeding  along  the  same  lines,  the  next  set  of  four 
trial  rays  shrink  the  impact  area  to  typically  a  5  km  square.  The 
'homing  in*  ray  is  determined,  finally,  on  the  basis  of  the  last  set 
of  four  impact  points.  The  ray  is  considered  to  have  successfully 
'homed  in'  if  it  returns  to  within  1  km  of  the  receiver  location. 

The  'homing  in'  ray  for  an  initial  time  serves  as  a  guide  to  arrive 
at  the  first  guess  trial  rays  for  a  subsequent  time.  Figure  4-1 
illustrates  the  sequence  of  trials  used  to  accomplish  the  homing  of 
an  extraordinary  ray  of  16.078  MHz  transmitted  at  70  seconds  after 
the  release  of  the  cloud.  It  was  found  that,  in  practice,  it  is  often 
easier  to  guess  the  bound  rays  in  elevation  than  in  azimuth  and  con¬ 
sequently  instead  of  (3x3),  a  (2  x  3)  initial  matrix  of  rays  as 
shown  in  the  figure  was  found  sufficient  to  enclose  the  receiver.  In 
selecting  the  next  set  of  four  rays,  some  account  is  made  for  the 
skewness  of  the  elevation  and  azimuth  lines  instead  of  choosing  the 
angles  associated  with  four  evenly  spaced  points  around  the  receiver. 
The  final  set  of  four  trial  rays,  however,  are  selected  evenly  to 
center  the  receiver.  The  sequence  shown  in  the  figure  represent  a 
shrink  in  the  elevation  angle  increment  from  0.08  to  0.0008°  and 
azimuthal  increment  from  0.3  to  0.0052.  On  obtaining  the  homing  ray, 
the  scattering  cross  section  presented  by  the  cloud  to  this  ray  is 
determined  as  described  below. 

4.1.1  RADAR  CROSS  SECTION  OF  THE  CLOUD 

The  scattering  cross  section  of  the  cloud  <7,  defined  in  terms 
of  an  equivalent  isotropic  scatter,  is  given  by  the  bistatic  radar 
equation  as: 
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Figure  4-1  Illustrates  the  Method  for  Determination  of  the  'Homing' 
Ray,  A,B,  C,D,  etc.  are  the  Ground  Impact  Points  When  a 
Bundle  of  Rays  Incremented  in  Azimuth  and  Elevation  are 
Transmitted 


a  = 


16 tt2  R  V2 

P  G4-  A 
t  t  r 


(4-1) 


where 


G,  = 


R,  « 


R  = 


'  A  = 


Transmitted  power 

Effective  gain  of  the  transmitting  antenna 
Range  to  the  scattering  volume  from  the  transmitter 
Range  to  the  scattering  volume  from  the  receiver 
Received  power 

Effective  collecting  area  of  the  receiving  antenna 


2  .  . 

A  can  be  expressed  as  Gra  / 4v ,  where  Gr  is  the  effective  gam  of  the 

receiving  antenna  and  A  is  the  wavelength  of  the  received  signal.  The 
scattering  cross  section  of  the  cloud  can  be  calculated  from  Equation 
4-1  following  the  technique  developed  by  Croft  (1967)  and  generalized 
recently  by  Georges  and  Stephenson  (1969)  for  3D  ray  tracing  simula¬ 
tion  of  hf  ground  backscatter  signatures.  The  method  derives  an 
expression  for  the  received  power  P^  which  results  in  an  expression 
for  a  free  of  system  parameters.  This  makes  it  possible  to  calculate 
a  solely  on  the  basis  of  the  ray  tracing  results.  The  method  in 
brief  is  as  follows: 


A  bundle  of  rays  are  transmitted  about  the  direction  of  the 
homing  ray  in  a  small  solid  angle  defined  by  an  elevation  angle 
increment  of  A/3  and  azimuth  angle  increment  of  A4> .  The  flux  of  energy 
leaving  in  this  solid  angle  will  refract  through  the  cloud  and  strike 
the  ground  on  return  over  an  area  denoted  as  A.  The  amount  of  power 
reaching  the  impact  area  is  a  fraction  K  times  the  power  fed  into  the 
flux  tube  of  transmission  where  the  fraction  is  determined  by  the 
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dissipation  losses  in  the  medium.  When  the  transmission  is  centered 
at  an  elevation  angle  of  then  the  solid  angle  into  which  the  energy 
is  injected  becomes  A|3  A1!5  cos  J3.  The  power  leaving  this  so_id  angle 

is  given  as: 

P^  G  A|3  A<i>  cos  P. 

t _ £ - - -  (4-2) 


when  this  power  approaches  the  impact  area  A  at  the  elevation  angle 
then  the  received  power  is  given  as: 


P 

r 


P.  G.  Af3A^  cos  p  A  K 

_fc _ £ - - - - - i - 

4tt  A  sin 


(4-3) 


The  fraction  K  is  given  as  exp  (-X/4.343)  where  X  is  the  total 
absorption  in  dB  suffered  by  the  homing  ray.  The  Equation  4-3 
for  P^  applies  strictly  to  a  cw  radar  and  in  the  case  of  a  pulse  radar 
it  should  be  multiplied  by  a  correction  factor  to  take  account  of  the 
pulse  spreading  in  time.  If  T  and  t  +  At  are  the  pulse  lengths  of  the 
transmitted  and  received  signals  the  correction  factor  is  given 
approximately  as  (t/t+At) .  This  correction  is  considered  small  and 
neglected. 

Substituting  for  P^  in  Equation  4-1  leads  to 


7 


2  2 

4r  R1  R2 


A(5  A4>  cos  (3  K 


A  sin  i 


(4-4) 


In  the  Equation  4-4  for  ■r,  the  parameters  P,  A/3  and  A4>  are  the  specified 
input  to  the  ray  tracing  program.  The  output  provides  K  and  all  the 
information  necessary  to  compute  R^,  R^  and  A.  The  geometry  involved 
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to  calculate  A  is  the  same  as  described  by  Goerges  and  Stephenson 
(1969)  and  is  as  shown  in  Figure  4-2.  Consider  a  set  of  four  rays 
that  constitute  a  rectangular  flux  tube.  The  rays  while  passing 
through  the  cloud  undergo  lateral  deviations  and  on  hitting  the 
ground  form  a  spherical  quadrilateral.  The  area  of  the  quadrilateral 
can  be  calculated  from  the  coordinates  of  the  four  corners  supplied  by 
the  ray  tracing  program.  In  practice,  the  four  rays  of  the  flux  tube 


Figure  4-2  Geometry  of  an  Electromagnetic  Flux  Tube  Illuminating  a 
Patch  of  Ground  on  Reflection  from  the  Plasma  Cloud 


are  so  closely  packed  that  the  quadrilateral  can  be  considered  plane 
and  a  simple  formula  for  the  area  can  be  applied.  If  r  ,  r  ,  r  and 

A  “  J 

r^  are  the  ground  ranges  to  the  four  landing  points  and  6  ,  9^,  9 „ 


and  0.  are  their  azimuths  then. 

4 

n  _  _L_  r3  ~  rl  r4  "  r: 

*  sin  a  sin  7 


sin  («  +  7) 


(4-5a) 


where 


tan  a  = 


tan  7 
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(4-5b) 


(4-5c) 


and  is  the  radius  of  the  Earth.  This  amounts  to  calculating  the 
area  of  a  quadrilateral  by  one-half  the  product  of  its  diagonals 
times  the  sine  of  the  angle  between  them.  The  area  is  calculated 
using  the  final  set  of  trial  rays  such  as  shown  in  Figure  4-1  provided 
it  is  less  than  about  5  km  square  and  the  receiver  is  reasonably  well 
centered.  Otherwise,  a  new  set  of  four  rays  closely  spaced  about 
the  homing  ray  are  used  to  meet  the  requirement. 

The  ranges  and  R^  to  the  reflection  point  from  the  transmitter 
and  the  receiver  are  calculated  by  forming  two  triangles  as  shown 
in  Figure  4-3.  Considering  the  triangle  TSO,  the  range  R^  is 


obtained  as: 


Ri  - 


Rr  sin  cos^  +  (3) 


(4-6) 


The  ground  range  to  the  subreflection  point  TS '  taken  from  the  ray 
tracing  output  provides  the  central  Earth  angle  £  . 


The  range  R^  is  obtained,  considering  the  Triangle  RSO,  as: 
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R2  "  K  sin  y 2  +  !RE  ■  Re  sin  e2  +  h) 2 


(4-7) 


The  ground  range  RSZ  which  provides  %2  is  calculated  from  the 

spherical  triangle  TRS^  knowing  the  ranges  TS  and  ^  and  the  azimuth 
of  the  subreflection  point  with  respect  to  TR. 

The  scattering  cross  section  cr  is  calculated  by  substituting  the 
Equations  4-5,  4-6  and  4-7  in  4-4. 


reflection  point 


0  CENTER  OF  EARTH 

Figure  4-3  Illustration  to  Evaluate  the  Ranges  R  and  »  to  the 

Reflection  Point  from  the  Transmitter  and  the  Receiver 
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4.2  HARD  EXPANDING  ELLIPSOID  METHOD 
4.2.1  RADAR  CROSS  SECTION  (RCS) 

The  radar  cross  section  of  a  smooth  body  is  given  from  considera¬ 
tions  of  geometric  optics  as  (Crispin  and  Maffett,  1968) : 


a  ~  v  rir2 


(4-8) 


Where  R^  and  R2  are  the  principal  radii  of  curvature  at  the  specular 
reflection  point  and  are  both  assumed  to  be  large  in  comparison  with 
the  operating  wavelength.  If  the  equation  of  the  target  surface  is 
expressed  in  the  form  z  =  f(x,y),  then 


R1R2  = 


(1  +  f2  +  f  2)2 
_ X _  V 


f  f  f 
xx  xy  -  xy 


(4-9) 


where  f  =  f  -  f 

x  3x  y  3y  xx 
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The  surface  of  interest  here  is  that  of  an  ellipsoid  for  which  the 
equation  can  be  written  as: 


(4-10) 


Using  the  coordinate  system  shown  in  Figure  4-4,  the  following  expres¬ 
sion  can  be  derived  for  monostatic  radar  cross  section  of  the  general 
ellipsoid. 


_ 7T  a  b  c _ 

J  "  [  2  .  2  ,  22  .  2  .2  2  2  1 2  (4-11) 

[a  sm  6  cos  0  +  b  sm  9  sin  <p  +  c  cos  6 J 

The  Equation  4-11  for  the  case  of  Barium  ion  clouds,  which  can  be 
approximated  to  a  prolate  spheroid  (a  =  b) ,  can  be  reduced  to 
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Figure  4-4  Coordinate  System  Adopted  to  Obtain  An  Expression  for 
the  Radar  Cross  Section  of  the  Ellipsoid 


a  = 


.4  2 
7 r  b  c 


b2  sin2  0  +  C2  cos2  e}2 


(4-12) 


The  radar  configuration  in  the  experiment  of  interest  here  (test 
Mulberry)  is  approximately  monostatic  and,  therefore,  the  Equation  4-12 
can  be  used  to  calculate  the  cross  section.  The  parameters  c  and  b, 
the  major  and  minor  axes  of  the  cloud,  and  the  aspect  angle  p  to  the 
specular  point  from  the  radar  are  evaluated  at  any  given  time  as 
described  in  the  following  paragraphs. 


4. 2. 1.1  Determination  of  c  and  b 


The  electron  density  in  the  cloud  is  given  according  to  the 
Equation  3-2.  The  parameters  c  and  b  are  the  major  and  minor  axes 

ir  with  the  electron  densitv  critical  to  the 


i  _  . 


i»M#] 


d 


[•111 


operating  frequency  of  the  radar.  If  n  denotes  the  critical  density, 

P 

it  is  given  by  considering  the  variation  along  the  axis  (r  =0)  as: 


n  = 
P 


Nj.  (1  -  e  t//T)  exp.  [-c2/(h2  +  ^D^t)] 
tt3/2  (h2  +  4Dnt)1/2  (h2  +  4Dit) 


(4-13) 


From  Equation  4-13  one  obtains 


«  =  <ho  +  4Dllt)  1/2  ^  — vr4-— Sr-5 - 

n  tt '  (h  +  4D.,t)  ^  ^  (h  +  4D,t) 
L  p  O  H  o  1  - 

(4-14) 


e't/T) 


similarly  considering  the  variation  along  a  radius  vector  at  z  =  0, 
one  aets 


b  =  (h2  +  4D1t)1//2  In 


N.  -t/T, 

i  (1  -  e 


n  7r3/2  (h2  +  4D„t)1/2  (h2  +  4D.t) 
D'li 


(4-15) 


4. 2. 1.2  Determination  of  the  Aspect  Angle  9 


Consider  a  cloud  centered  Cartesian  system  (x,  y,  z)  as  defined 
in  Section  3  with  z  axis  anti-parallel  to  the  direction  of  the  Earth' s 
magnetic  field.  Let  the  coordinates  of  the  radar  in  this  system  be 
(x^,  y  ,  z^)  and  those  of  the  specular  point  be  denoted  as  (x^,  yQ, 
z  )  .  The  coordinates  of  the  specular  point  are  evaluated  by  con¬ 
sidering  a  unit  vector  normal  to  the  surface  at  that  point.  The 
equations  for  the  ellipsoid  and  the  unit  normal  vector  are: 


T  2  2  ( 

=  c  -  c  ( 


2  .  2  V 


(4-16) 


N  = 


if  +  j  f  -  k 
xo  yo _ 

(f2  +  f2  + 1)  1/2 

xo  yo 


(4-17) 
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are 


Where  i,  j,  k  are  the  unit  vectors  along  x,  y,  z  and  f  ,  f  , 

Qz  0z  xo  yo 

the  partial  derivatives  —  and  at  the  point  (x  ,  y  ,  z  )  .  Since 

8x  Sy  o  Jo  o 

the  normal  vector  is  along  the  line  joining  the  specular  point  to  the 
radar,  the  unit  vector  N  can  also  be  expressed,  using  the  geometry 
shown  in  Figure  4-5,  as: 


N- 

PR 


i  (xx  -  xo)  +  j  (yx  .  yQ)  +  k  {z1  -  o 


_{X1  “  Xo)2  +  {yl  "  yo)2  +  (zi  "  Zo)2_ 


2  1 1/2 


Equating  4-17  and  4-18  and  substituting  (~ 
2 

c 


2 

c  x 

— r  “)  for  f  and 
.2  xo 

b  z 

o 


(4-18) 


(-  —  — )  for  f  leads  to 

b2  z  y° 

o 


2 

c  z. 


X,  =  X  1  -  —  +  T 


y  =  1  - 


b2  z 

c 

2 

b2  z 


(4-19) 


(4-20) 


Figure  4-5  The  Geometry  to  Express  Unit  Vector  Normal  to  the  Cloud 
Surface,  at  the  Specular  Point  from  the  Radar 
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Equations  4-19  and  4-20  are  squared  and  then  added  to  yield,  on 
substitution  of  [b2  -  (b2/c2)z2]  for  x2  +  y2,  a  quartic  in  zq  of  the 


form 


Az4  +  Bz3  +  Cz2  +  Dz  +  E  =  0 
o  o  o  o 


(4-21) 


where 


',2  2 
b  c 

2  +  ,2  "  2 
•c  b 


B=  2^-2-^ 


r  2  4 

„  2  2  c  2  ,2  c^_2 

C  =  I xn  +  y.,  +  “  z1  -  b  -  ~  +  2c 


c  2 

D  =  2  ~  2  -  2c  z 

b  1  X 


2  21 


Since  the  dimensions  of  the  cloud  are  very  small  compared  to  the  dis¬ 
tance  between  the  radar  and  the  cloud,  a  good  approximation  to  the 
solution  of  the  Equation  4-21  can  be  obtained  by  dropping  the  higher 
order  terms  with  coefficients  A  and  B.  The  resulting  equation  is 

solved  for  z  and  of  the  two  oossible  values,  the  one  with  favorable 

o  ‘  ' 

geometry  for  specular  reflection  is  taken  and  substituted,  in  turn,  in 
Equations  4-19  and  4-20  to  get  xq  and  yQ.  The  coordinates  of  the  radar 
that  are  used  to  calculate  the  specular  point  are  given  by: 


,  *1  fair  cos  <j>  cos  y  sin  $  cos  y 
1  mo  Amo  mo  Ai 


yl  =  b 


‘sin^o  Xml 


-sin  4> 


cos  4> 


0  Y„ 


cos  4>  sin  y  sin  <j>  sin  y  cos  y 


mo  ml 


(4-22) 

The  elements  a,  b*  and  c*  and  the  transformation  matrix  are  the  same 

as  in  Equation  3-7.  (x  , ,  y  .  z  )  are  the  coordinates  of  the  radar 

ml  ml  ml 

in  the  geomagnetic  dipole  coordinate  system  defined  in  paragraph  3,2. 
Knowing  the  coordinates  of  the  specular  point,  the  aspect  angle  0  can 
be  calculated  by  taking  the  scalar  product  of  the  unit  vector  N  and  a 
unit  vector  along  the  magnetic  field  B. 


N  =  ^  i  +  £2  j  +  J3  £ 

B  =  mx  i  +  m2  j  +  m3  £ 


then 


cos  6  =  ^.m.  +  i_mn  +  iLnu 

1  1  2  2  3  3 

The  direction  cosines  are  given  as:  m,  =m  =0  m  =  1 

r  2  2  v1/2  2  '  3 

and  £  =  -  f  +  f  ■ 1  +  jl 

3  L  xo  yo  J 


* 

Not  to  be  confused  with  major  and  minors  axes  of  the  ellipsoid. 
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substituting  for  f  and 

xo 


f  one  obtains: 
yo 


cos  9 


/  2  2. 
(xn  *  yj 


+  1 


-1/2 


4.2.2  DETERMINATION  OF  DOPPLER  SHIFT 


(4-23) 


Let  (xol,  yQl,  zQl)  and  (xq2#  yo2/  z^)  be  the  coordinates  of 

the  specular  reflection  points  on  the  cloud  at  times  t  and  t  +  At 
as  shown  in  Figure  4-6.  The  corresponding  ranges  to  these  reflection 
points  from  the  radar  are  given  as: 


Figure  4-6  Illustrate  the  Position  of  the  Specular  Point  at  Two 

Instants  of  Time  to  Calculate  the  Path  Difference  to  it 
from  the  Radar 
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(4-24) 


Ri  =L(xi"x°i}  +  (Yi“yoi)  +  ‘vv  J 

Ez  =[(xrxo2)2  +  (yry02>2  +  urzo2)2]' 


(4-25) 


Since  the  signal  traverses  in  free  space  the  difference  between  the 
two  ranges  AR(=R2“R^)  represents  the  difference  in  the  phase  path 
length  to  the  reflection  point*  Hence  the  Doppler  shift  of  the 
reflected  signal  at  the  radar  is  obtained  by 


Af  = 


-2AR 

Mt 


(4-26) 


The  Doppler  shift  given  in  the  Equation  4-26  is  smeared  over  the 
time  interval  At.  The  time  increment  should,  therefore,  be  made  as 
small  as  possible  without  running  into  computer  round  off  errors  in 
calculating  AR.  Figure  4-7  shows  a  typical  case  to  illustrate  how 
the  Doppler  shift  approaches  its  true  value  in  an  asymptotic  manner 
as  At  is  decreased  from  0.5  to  0.01  seconds.  If  At  is  allowed  to 
take  lower  values,  it  would  result  in  a  fluctuation  in  Af.  Since 
the  lower  limit  of  At  is  time  dependent,  it  is  allowed  to  vary  over 
an  interval  at  each  time  and  the  asymptotic  value  such  as  shown  in 
the  Figure  4-7  is  picked  for  the  Doppler  shift. 

4.3  RESULTS  AND  DISCUSSION 


4.3.1  RAY  TRACING  RESULTS 

Figure  4-8  shows  ray  path  projection  plots  in  both  vertical  and 
ground  planes  at  100  seconds  after  release  for  all  the  three  modes  name¬ 
ly,  the  no-field  the  ordinary  and  the  extraordinary.  The  ray  path  seg¬ 
ments  close  to  the  reflection  level  where  there  is  an  appreciable  bending 
are  plotted  on  an  enlarged  scale  for  the  vertical  plane  in  Figure  4-9 
to  show  how  the  three  modes  refract  in  the  denser  regions  of  the  cloud. 


MULBERRY  ELLIPSOID  MODEL 


(  *H)  jy 
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Figure  4-7  Illustrate  the  Limiting  Process  Involved  in  Obtaining  the  Doppler 
Shift  from  the  Path  Difference  Calculation 


olicnn^o1!^  in  Both  Vertical  and  Ground 
0  Seconds  After  Release  for  the  Three  Modes 

dmary  and  Extraordinary,  at  F  =  10.31  mL  ' 
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MULBERRY.  T= 100.  F=  10-310 


Figure  4-9  The  Vertical  Ray  Path  Segments  Close  to  the  Reflection 
Level  to  Show  How  the  Three  Modes  Refract  Through  the 
Denser  Regions  of  the  Cloud 
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The  movement  of  the  reflection  point  with  respect  to  the  center  of 
the  cloud  and  the  refraction  suffered  by  the  rays  as  function  of  time 
are  shown  in  Figures  4 -10a,  b  and  c  and  4-lla  and  b  in  terms  of 
the  azimuth  and  elevation  angles  of  the  cloud  center,  the  reflection 
point  and  of  transmission.  It  can  be  seen  that  the  direction  of  the 
reflection  point  walks  away  from  the  center  during  the  expansion 
phase  as  the  cloud  becomes  elongated  and  walks  back  at  later  times 
as  the  reflection  surface  recedes  toward  the  center.  The  deviations 
in  azimuth  and  elevation  due  to  refraction  are  of  some  importance  and 
are  therefore  plotted  separately  in  Figure  4-12.  The  refraction  in  the 
case  of  no-field  ray  is  entirely  due  to  the  density  gradients  in  the 
ionization  at  levels  below  the  height  of  reflection.  The  deviation 
is  nearly  zero  at  the  beginning  and  increases  continuously  at  a  rate 
that  builds  up  rapidly  with  time.  The  0  and  X  rays  undergo  refrac¬ 
tion  that  is  partly  due  to  density  gradients  and  partly  due  to  the 
magnetic  field.  The  refraction  effects  due  to  density  gradients  and 
the  magnetic  field  are  found  to  act  in  the  same  direction  for  the 
X  ray  while  they  get  cancelled  partially  for  the  0  ray.  The  curves 
for  the  0  ray  pass  through  zero  indicating  the  time  when  the  two 
effects  cancel  each  other  exactly.  The  refraction  is  found  to  be  the 
highest  for  the  X  ray  with  azimuth  and  elevation  deviations  of 
1.65  and  -0.35  respectively  (the  A<j>  is  in  the  ground  plane  and 
should  be  multiplied  with  cos  (3  to  obtain  it  in  the  transmission 
plane).  The  cross  over  of  the  curves  for  the  no-field  and  the  X  ray 
is  due  to  the  fact  that  the  deviations  increase  most  rapidly  at  the 
time  of  penetration  which  comes  slightly  earlier  for  the  no-field  ray. 

The  defocussing  nature  of  the  cloud  is  illustrated  in  Figures 
4-13  to  4-16.  The  ray  plots  shown  are  for  the  frequency  of  10.31  MHz 
and  for  the  time  of  180  seconds  after  release.  The  plots  are  pro¬ 
jections  of  the  ray  paths  on  the  ground  and  the  vertical  plane 
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Figure  4-12  The  Azimuth  and  Elevation  Angle  Deviations  Suffered 
by  the  Three  Modes  Due  to  Refraction  up  to  the 
Reflection  Level  at  F  =  10.31  MHz 
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MULBERRY. Trl8Q.FrlQ. 310 
XTRRORDINRRY 

ELEVATION  SCAN  -  VERTICAL  PLANE 


Figure  4-13  Ray  Path  Projections  in  Vertical  Plane  of  a  Fine  Bundle 
of  Eleven  Rays  (Appears  as  Thick  Line  on  the  Plot) 
Transmitted  Over  an  Elevation  Scan  Interval  of  0.02  . 
The  Spread  of  the  Rays  on  Reflection  Illustrate  the 
Defocussing  Effects  of  the  Cloud 
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Figure  4-14  Ray  Path  Projections  in  Vertical  Plane  of  a  Fine  Bundle 
of  Eleven  Rays  (Appears  as  Thick  Line  on  the  Plot) 
Transmitted  over  an  Azimuth  Scan  of  Size  0.025  .  The 
Spread  of  the  Rays  on  Reflection  Illustrate  the 
Defocussing  Effects  of  the  Cloud 
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MULBERRY .T=180  ,F  =  10  .310 
XTRRORD  INflRY 

ELEVATION  SCAN-HORIZONTAL  PLANE 


th  Projections  on  the  Ground  Plane  of  a  Fine  Bundle  of  Eleven  Rays 
Line  on  the  Plot  from  TX)  Scanned  in  Elevation  by  0.02  Degree 


Figure  4-16  Ray  Path  Projections  on  the  Ground  Plane  of  a  Fine  Bundle 
of  Eleven  Rays  (Thick  Line  on  the  Plot  from  TX)  for  an 
Azimuth  Scan  of  Size  0.025° 


51 


connecting  the  transmitter  and  the  receiver.  They  describe  how  a 

narrow  beam  of  incident  energy  confined  either  to  an  elevation  plane 

or  an  azimuth  plane  become  defocussed  on  reflection  from  the  cloud. 

The  beam  sizes  A4>  cos  0  and  A f3  represented  by  the  scans  in  azimuth 

and  elevation  are  respectively  0.025  degrees  and  0.02  degrees.  They 

are  so  narrow  in  extent  that  they  could  be  seen  only  as  a  thick  line 

on  the  plots.  The  energy  dispersion  on  the  ground  is  found  to  be  to 

an  extent  of  about  130  km  and  160  km  corresponding  to  the  elevation 

and  azimuth  scans  of  the  ground  plane  plots.  The  spreading  of  the 

rays  due  to  the  defocussing  effects  of  the  cloud  in  this  case  are 

5 

greater  by  a  factor  of  about  8  x  10  relative  to  free  space  spreading. 
The  defocussing  effects  are  not  uniform  over  the  entire  illuminated 
area  since  the  rays,  while  equally  incremented  in  the  scan,  seem  to 
bunch  closer  in  certain  areas.  It  should  be  mentioned  that  some  of 
the  bunching  effect  seen  especially  with  the  azimuthal  scan  in  the 
vertical  plane  is  due  to  the  distortion  resulting  from  the  projection. 
The  cloud  is  found  to  become  increasingly  more  defocussing  with  time 
as  it  gets  more  and  more  elongated.  On  the  basis  of  the  impact  area 
calculation  described  in  paragraph  4.1,  it  is  found  that  the  defocuss¬ 
ing  factor  (relative  to  free  space)  for  the  X  ray  at  10.31  MHz 
.  3  7 

increases  from  7.6  x  10  to  1.85  x  10  over  the  interval  5  to  190 

4  8 

seconds.  It  will  increase  from  2.27  x  10  *  to  1.47  x  10  in  the 
interval  10  to  85  seconds  for  16.078  MHz. 

The  Doppler  and  the  scattering  cross  section  results  synthesized 
by  the  ray  tracing  technique  are  presented  in  Figures  4-17  and  4-18 
for  the  two  operating  frequencies  10.31  MHz  and  16.078  MHz.  In 
addition  to  the  ordinary  and  extraordinary  modes,  results  have  been 
obtained  also  for  the  no-field  case  for  10.31  MHz  to  study  the 
magnetic  field  effects.  The  shapes  of  the  curves  showing  the  variation 
of  the  Doppler  and  the  cross  section  with  time  are  found  to  be 
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Figure  4-17  The  Time  History  of  the  Doppler  Shift  and  the  Scattering 
Cross  Section  Synthesized  by  the  Ray  Tracing  Method  for 
No~Field,  Ordinary  and  Extraordinary  Modes  at 
F  =  10.31  MHz 
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Figure  4-18  The  Time  History  of  the  Doppler  Shift  and  the  Scattering 
Gross  Section  Synthesized  by  the  Ray  Tracing  Method  for 
Ordinary  and  Extraordinary  Modes  at  F  =  16.078  MHz 


identical  lor  all  the  three  modes.  The  Doppler  and  the  ctc<ls. 

variations  for  the  X  mode  lag  in  time  to  that  of  the  0  mode.  The  lag 

in  Doppler  widens  as  the  cutoff  times  are  approached.  The  variations 

for  the  ‘no  field*  (N)  fall  in  between  that  of  the  0  and  X,  but  closer 

to  the  0  mode.  The  differences  between  the  0  and  X  modes  are  primarily 

due  to  the  fact  that  they  are  supported  by  two  plasma  frequency 

surfaces  whose  expansion  rates  are  changing  relative  to  each  other 

with  time.  The  plasma  frequency  contour  supporting  the  0  and  N  modes 

being  the  same,  the  differences  between  these  two  indicate  the  effects 

due  to  the  magnetic  field.  At  the  operating  frequency  of  10.31  MHz 

the  Doppler  shift  falls  from  17  to  -5  Hz  in  the  time  interval  of 

5  seconds  after  release  to  the  cutoff.  The  cross-over  from  positive 

to  the  negative  Doppler  occurs  roughly  halfway  in  the  time  interval, 

7  2 

The  scattering  cross  section  reaches  a  maximum  of  about  2  x  10  m 
within  20  to  30  seconds  after  release  and  thereafter  falls  continuously 
at  a  rate  that  increases  rapidly  with  time.  It  implies  that  the 
signal  drops  out  rather  abruptly.  The  results  for  16.078  MHz  share 
all  the  basic  features  observed  for  10.31  MHz.  The  peak  scattering 
cross  section  at  16.078  MHz  is  found  to  be  smaller  than  at  10.31  MHz 
by  a  factor  5  for  the  ordinary  ray  and  4  for  the  extraordinary  ray, 

4.3.2  COMPARISON  WITH  ELLIPSOID  MODEL 

The  Doppler  and  the  scattering  cross  section  have  been  calculated 
also  for  the  ,'ield  aligned  ellipsoid  model  using  the  method  described 
in  paragraph  4.2.  The  results  are  presented  in  Figures  4-19  and 
4-20.  The  curves  denoted  by  0  and  X  in  these  figures  represent  two 
slightly  different  operating  frequencies  which  correspond  to  the 
plasma  frequencies  appropriate  to  the  ordinary  and  the  extraordinary 
modes.  The  comparison  between  the  ellipsoid  model  and  the  ray  tracing 
is  intended  to  bring  out  the  effects  on  the  Doppler  and  the  scattering 
cross  section  due  to  magnetic  field  and  the  ionization  lying  below 
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Figure  4-19 
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The  Time  History  of  the  Doppler  Shift  and  the  Scattering 
Cross  Section  Synthesized  by  the  Hard  Ellipsoid  Method 
for  F  =  10.31  and  9.902  MHz.  (The  Two  Frequencies 
Correspond  to  the  Plasma  Frequencies  that  Support  the 
0  and  X  Modes) 
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Figure  4-20  The  Time  History  of  the  Doppler  Shift  and  the  Scattering 
Cross  Section  Synthesized  by  the  Hard  Ellipsoid  Method 
for  F  =  16.078  and  15.673  MHz.  (The  Two  Frequencies 
Correspond  to  the  Plasma  Frequencies  that  Support  the 
0  and  X  Modes) 


the  reflection  level.  The  general  shapes  of  the  curves  showing  the 
temporal  variations  of  the  two  parameters  are  similar  for  both  the 
models.  Considering  in  detail,  the  Doppler  variations  are  signifi¬ 
cantly  steeper  for  the  ellipsoid  model  than  for  the  ray  tracing.  The 
negative  phase  of  the  Doppler  is  much  more  pronounced  for  the  ellip¬ 
soid  model  with  the  transition  from  positive  to  negative  occurring 
significantly  earlier  than  for  ray  tracing.  The  transition  occurs 
for  the  ellipsoid  model  at  about  the  same  time  as  the  scattering 
cross  section  begins  to  decrease  whereas  it  lags  a  short  while  with 
ray  tracing.  This  suggests  that  an  effect  of  the  underlying  ioniza¬ 
tion  is  to  hold  the  Doppler  positive  for  a  short  duration  even  after 
the  contraction  phase  of  the  reflection  surface  is  initiated.  It 
causes  also  the  tailing  effect  of  the  Doppler  to  be  significantly 
inhibited.  The  defocussing  effects  of  the  underlying  ionization  which 
enhance  with  time  cause  a  remarkable  reduction  in  the  scattering  cross 
section  of  the  cloud.  The  ray  tracing  technique  yields  a  cross  section 
which  is  smaller  by  a  factor  that  increases  with  time  from  3  at  the 
peak  to  30  at  the  time  just  prior  to  penetration.  The  comparison 
thus  leads  to  the  conclusion  that  the  simple  hard  ellipsoid  model  is 
inadequate  for  an  accurate  synthesis  of  the  hf  radar  signatures. 
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SECTION  5.  EXPERIMENTAL  RESULTS 


5.1  TEST  GEOMETRY  AND  INSTRUMENTATION 

The  Pre-SECEDE  Series  of  Barium  release  tests  were  conducted 
during  the  fall  of  1969  at  White  Sands  Missile  Range,  New  Mexico.  T 
releases  were  observed  by  nine  hf  radar  systems  of  which  eight 
(6  pulse  and  2  cw)  had  been  operated  by  Raytheon  with  transmitters 
located  at  Paxton  Siding  and  receivers  located  at  Twin  Butes. 

Figure  5-1  shows  the  geometry  of  the  radar  sites  with  respect  to  the 


w 

Figure  5-1  Illustration  to  Show  the  Cloud  Geometry  With  Respect 
to  the  Transmitting  (TX)  and  Receiving  (RX)  Stations 


position  of  the  plasma  cloud  at  release  for  the  test  Mulberry.  The 
discussion  of  the  results  presented  in  this  section  is  limited  to 
the  pulse  measurements  made  at  10.31  and  16.078  MHz  for  which  the  ray 
tracing  synthesis  has  been  done. 

The  radars  employed  in  the  experiments  are  the  conventional  Phase 
Path  Sounders  with  the  transmitting  and  the  receiving  systems  phase 
locked  to  stable  frequency  standards  at  each  site.  This  enables  both 
phase  and  amplitude  of  the  received  signal  to  be  obtained  from  the 
observations.  The  pulse  radars  had  been  operated  at  a  40  pps  rate 
with  a  pulse  width  of  150  |is.  For  receiving,  a  bandwidth  10  kHz 
was  used  and  the  output  signal  was  recorded  on  analog  magnetic  tape 
at  a  15  kHz  i-f  level  for  later  processing.  The  antennas  used  for 
transmitting  10.31  and  16.078  MHz  were  vertically  directed  deltas. 

At  the  receiving  site,  a  delta  antenna  was  used  for  receiving  at 
10.31  MHz  and  a  vertically  directed  log  periodic  antenna  for 
receiving  at  16.078  MHz. 

5.2  EXPERIMENTAL  DATA 

The  experimental  data  recorded  on  the  tape  is  reduced  to  the  form 
of  doppler  spectra.  The  signature  display  format  is  such  that  the 
time  runs  from  left  to  right,  Doppler  shift  runs  from  bottom  to  top 
and  the  strength  of  the  echo  is  represented  by  the  darkness  of  the 
trace.  The  signal  processing  involves  phase  detection  of  the  i-f 
signal  run  off  from  the  tape  which  is  then  range  gated  and  sampled 
at  the  pulse  repetition  rate.  The  sampled  signal  results  in  a  box- 
carved  waveform  which  is  fed  into  a  UA-6A  spectrum  analyzer  to  obtain 
the  desired  Doppler  signature  display.  The  unambiguous  window  of  the 
Doppler  spectrum  is  equal  to  half  the  pulse  repetition  frequency  which 


in  this  case  is  20  Hz,  The  Doppler  components  that  fall  outside  the 
range  fold  back  into  the  window  and  cause  images  to  appear  along  with 
the  real  spectral  components.  Figure  5-2  shows  the  Doppler  signatures 
obtained  on  10.31  and  16.078  MHz  for  the  test  Mulberry.  The  true 
signature  is  simply  a  single  monotonically  descending  tone  and  during 
the  first  few  seconds  when  the  Doppler  falls  outside  the  unambiguous 
window  an  image  appears  as  labeled  on  the  figure.  The  frequency 
folding  takes  place  about  10  Hz  instead  of  20  Hz  since  the  signal  is 
phase  detected  with  an  offset  frequency  of  10  Hz.  The  periodic  notches 
that  appear  on  the  traces  were  identified  as  due  to  Faraday  fading  of 
the  received  signal. 

The  spectrally  analyzed  signal  can  be  displayed  also  in  the  form 
of  amplitude  spectral  cuts  in  time  to  measure  the  received  power  of 
the  signal  and  thereby  determine  the  scattering  cross  section  of  the 
cloud  as  function  of  time.  The  sequence  of  spectra  obtained  have  a 
time  resolution  of  1.5  seconds  for  both  10.31  and  16.078  MHz.  Figures 
5-3  and  5-4  show  the  spectral  sweeps  for  selected  times  for  the  two 
frequencies.  The  signal  spikes  appear  in  pairs  because  of  the  imaging 
of  the  spectral  components  about  the  folding  frequency  of  10  Hz.  The 
signal  amplitudes  are  subject  to  nonuniform  weighting  in  frequency 
imposed  by  the  method  of  sampling.  The  weighting  function  is 
sin  X/X  where  X  is  equal  to  o)T/2  while  co  is  the  angular  frequency  and 
t  is  the  sampling  time.  The  amplitudes  scaled  from  the  spectra  are 
corrected  for  the  nonuniform  response.  A  test  signal  injected  into 
the  receiver  at  the  operating  frequency  in  5  dB  steps  from  a  reference 
power  level  is  used  for  calibrating  the  amplitude  scale.  Figure  5-5 
shows  the  calibration  curves  constructed  for  10.31  and  16.078  MHz  for 
converting  the  arbitrary  amplitude  readings  to  absolute  signal  power. 
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Figure  5-3  Selected  Scans  of  the  Amplitude  Spectra  of  the  Signal  at 
F  =  10.31  MHz  for  Test  Mulberry.  (The  Folding  of  the 
Spectral  Components  Lying  Outside  the  Unambiguous  Doppler 
Window  Causes  Images  to  Appear  on  the  Display.) 
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Figure  5- 


4  Selected  Scans  of  the  Amplitude  Spectra  of  the  Signal  at 
F  =  16.078  MHz  for  Test  Mulberry.  (The  Folding  of  the 
Spectral  Components  Lying  Outside  the  Unambiguous  Doppler 
Window  Causes  Images  to  Appear  on  the  Display.) 
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Figure  5-5  The  Calibration  Curves  Used  to  Convert  the  Arbitrary  Units 
of  Amplitude  to  Absolute  Signal  Power 
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The  amplitude  and  power  plots  for  the  two  frequencies  are  presented 
in  Figures  5-6  and  5-7.  The  dashed  horizontal  lines  shown  on  the 
plots  correspond  to  the  average  noise  level  and  the  minimum  detectable 
signal  level  which  is  arbitrarily  set  at  one  standard  deviation  above 
the  noise  level.  The  systematic  fading  observed  in  the  signal  ampli¬ 
tude  is  considered  as  due  to  the  interference  between  the  ordinary  and 
extraordinary  modes  of  the  echo.  The  power  plots  are  constructed  from 
the  envelopes  of  the  amplitude  curves  and  hence  represent  the  sum  of 
the  powers  returned  in  the  two  modes. 


5.3  RESULTS 


The  power  measurements  of  the  received  signal  are  converted  to 
the  scattering  cross  sections  of  the  plasma  cloud  using  the  radar 
equation 
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The  parameters  of  Equation  5-1  are  the  same  as  defined  in  paragraph 

4.1.1.  The  gains  G  and  G  include  the  system  losses  and  are  given 

t  r 

as  (g  7  )  and  (g  7  )  where  g's  and  7*s  represent  respectively  the 
antenna  gains  and  the  efficiency  factors.  The  system  parameters  used 
to  calculate  the  scattering  cross  sections  are  listed  in  Table  5-1. 


TABLE  5-1 
SYSTEM  PARAMETERS 


A  (Meters) 

Pt  (Watts) 

gt 

gr 

7t 

7 

r 

18.66 

1.24  x  104 

4.60 

6.31 

0.295 

0.617 

29.40 

1.50  x  104 

3.70 

4.39 

0.490 

0.414 

66 


1  I  4(  |  I  I  I  I  I  I  I  I - T 


0143  20  [40  0146  20  40  014?  20  40  0t48  20  40  0149 

ft  TIME  (GMT! 

Figure  5-6  The  Amplitude  (Arbitrary  Scale)  and  the  Power  Variations 

of  the  Signal  at  10.31  MHz.  The  Power  Curve  is  Constructed 
from  the  Envelope  of  the  Amplitude  Plot.  The  Horizontal 
Dashed  Lines  Indicate  the  Average  Noise  and  the  Minimum 
Detectable  Signal  Levels. 
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Figure  5-7  The  Amplitude  (Arbitrary  Scale)  and  the  Power  Variations 
of  the  Signal  at  F  =  16.078  MHz.  The  Power  Curve  is  Con¬ 
structed  from  the  Envelope  of  the  Amplitude  Plot.  The 
Horizontal  Dashed  Lines  Indicate  the  Average  Noise  and  the 
Minimum  Detectable  Signal  Levels. 
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The  antenna  gains  in  the  direction  of  the  cloud  for  the  deltas  are 

obtained  by  computing  the  patterns  for  the  configurations  used  in  the 

experiments.  An  estimate  of  8  dB  is  adopted,  however,  for  the  log 

periodic  antenna.  The  efficiency  factors  7 ^  and  7  account  for  the 

t  r 

cable  and  impedance  mismatch  losses.  The  ranges  and  R2  to  the 
reflection  point  of  the  cloud  from  the  transmitter  and  the  receiver 
are  taken  from  the  results  of  the  ray  tracing  program.  Figure  5-8 
shows  how  the  ranges  vary  with  time  as  the  reflection  point  moves 
out  and  then  recedes  toward  the  center  of  the  cloud.  The  scattering 
cross  sections,  calculated  on  the  basis  of  the  above  parameters 
and  the  power  measurements  of  the  received  signals,  are  presented  for 
the  two  frequencies  in  Figures  5-9  and  5-10  along  with  the  Doppler 
information  read  out  of  the  amplitude  spectra.  The  minimum  detectable 
cross  section  a  is  also  indicated  on  the  plots  and  it  is  about 

3  X  10  rn  for  10.31  MHz  and  4.4  x  102  m2  for  16.078  MHz.  The  initial 
maximum  observed  in  the  scattering  cross  sections  on  the  two  frequen¬ 
cies  is  believed  to  be  due  to  thermal  ionization  that  appears  immedi¬ 
ately  on  release.  The  following  minimum  which  was  reported  to  be 
more  pronounced  at  higher  frequencies  is  caused  apparently  by  the 
decay  of  the  thermal  ionization  before  the  photoionization  builds 
up  significantly.  The  maximum  that  is  associated  with  the  photo - 
ionization  occurs  at  about  110  seconds  for  10.31  MHz  and  50  seconds 
for  16.078  MHz.  The  observed  Doppler  behavior  is  characteristic  of 
an  expanding  cloud  with  the  rate  of  expansion  decreasing  continuously 
with  time.  Further  discussion  of  the  results  is  deferred  to  the  next 
section  where  they  are  compared  with  those  synthesized  by  ray  tracing 
through  Gaussian  clouds. 
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Figure  5-10  Time  History  of  the  Doppler  Shift  and  the  Scattering 
Cross-Section  Derived  from  the  HF  Measurements  at 
F  =  16.078  MHz 
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SECTION  6.  COMPARISON  BETWEEN  EXPERIMENTAL  AND  RAY  TRACING  RESULTS 


The  time  dependent  Gaussian  model  adopted  for  the  ray  tracing, 
while  permitted  to  change  its  origin  from  one  time  to  another  on  the 
basis  of  the  optical  trajectory  data,  had  been  kept  stationary  at  any 
given  time  as  the  rays  were  traced.  Since  the  cloud  was  observed  to 
drift  continuously  in  position,  there  is  a  Doppler  component  associated 
with  this  bodily  motion  which  should  be  added  to  that  computed  for  the 
expanding  stationary  cloud  to  allow  a  comparison  between  the  theory 
and  the  observations.  A  computer  program  is  available  which  calculates 
the  drift  Doppler  when  the  coordinates  of  the  cloud  as  a  function  of 
time  and  the  operating  frequency  are  given.  Figure  6-1  shows  the 
drift  Dopplers  as  function  of  time  at  10.31  and  16.078  MHz  for  the 
test  Mulberry.  Using  these  values,  the  computed  Dopplers  have  been 


Figure  6-1  The  Doppler  Shift  Associated  With  the  Drift  of  the  Plasma 
Cloud  for  F  =  10.31  and  16.078  MHz 
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corrected  and  are  presented  along  with  the  scattering  cross  sections 
and  the  corresponding  experimental  results  in  Figures  6-2  and  6-3. 

The  synthesized  curves  of  the  extraordinary  mode  are  selected  for 
comparison  since  they  represent  the  actual  conditions  at  the  time  of 
the  signal  drop-out  which  is  of  maximum  importance.  There  is  an 
apparent  discrepancy  in  the  cutoff  times  indicated  by  the  experimental 
and  the  synthesized  curves.  This  seems  to  be  due  to  the  fact  that 
the  parameters,  the  total  ion  content  and  the  diffusion  coefficients, 
which  determine  the  lifetime  of  the  cloud  for  different  operating 
frequencies  have  been  kept  constant  in  ray  tracing  synthesis  where 
in  fact  they  change  with  time.  The  discrepancy  can  be  removed  by 
adjusting  the  time  scales  of  the  synthesized  curves  to  match  with  that 
of  the  experimental  curves.  The  adjustment  is  not  made,  however,  in 
Figures  6-2  and  6-3. 

The  Doppler  and  the  scattering  cross  section  synthesized  by  ray 

tracing  are  found  to  differ  quite  significantly  from  those  measured 

experimentally.  The  differences  are  consistent  for  both  the  operating 

frequencies,  10.31  and  16.078  MHz.  The  most  important  difference  in 

the  Doppler  behavior  is  that  the  synthesized  Doppler  based  on  a 

Gaussian  model  has  a  tail  at  the  time  of  signal  drop-out  which  is  not 

experimentally  observed.  It  has  been  speculated  that  the  signal  might 

be  below  the  level  of  system  sensitivity  and  could  have  been  lost  in 

the  background  noise  at  the  time  when  the  Doppler  tail  is  expected 

to  be  seen.  The  estimated  minimum  detectable  cross  section  a  .  and 

mxn 

the  ray  tracing  results,  however,  clearly  show  that  the  signal  level 

is  above  the  system  sensitivity  and  the  doppler  tail  should  have  been 

observed  if  the  distribution  were  truly  Gaussian.  Even  if  a  . 

mxn 

were  greater  by  two  orders  of  magnitude  than  estimated,  still  the  signal 
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Figure  6-2  Comparison  Between  The  Experimental  And  The  Synthesized 
Time  Histories  Of  The  Doppler  Shift  And  The  Scattering 
Cross  Section  At  F  =  10.31  MHz 
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MULBERRY  F  =  16.078  MHz 
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Figure  6-3  Comparison  Between  the  Experimental  And  The  Synthesized 
Time  Histories  Of  The  Doppler  Shift  And  The  Scattering 
Cross  Section  At  F  =  16.078  MHz 
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should  be  detectable  over  a  significant  portion  of  the  Doppler  tail. 
Since  the  uncertainties  in  the  system  losses  are  not  likely  to  be 
greater  than  20  dB,  the  ray  tx*acing  results  lead  to  the  conclusion 
that  the  electron  density  distribution  is  unlike  Gaussian  during  the 
late  time  evolution  of  the  cloud.  Considering  the  general  behavior, 
the  experimental  Doppler  is  found  to  be  consistently  greater  than  that 
synthesized.  It  was  noted  earlier  from  a  comparison  of  the  Dopplers 
by  hard  ellipsoid  and  ray  tracing  methods  that  the  effect  of  the 
'softness'  of  the  cloud  is  to  contribute  some  positive  Doppler.  In 
analogy  to  that,  the  comparison  between  the  experimental  and  the  ray 
tracing  Doppler  suggests  that  the  cloud  is  softer  than  that  implied 
by  the  adopted  Gaussian  distribution.  The  scattering  cross  sections 
synthesized  by  ray  tracing  are  10  to  20  dB  greater  than  those  derived 
from  the  experimental  data.  Since  the  uncertainties  in  antenna  gains 
and  system  losses  are  not  likely  to  add  up  to  this  magnitude,  the 
true  cross  sections  of  the  cloud  seem  to  be  lower  than  that  predicted 
by  the  Gaussian  model.  This  suggests  that  the  cloud  is  softer  than  the 
Gaussian  causing  a  greater  defocusing  and  thereby  a  reduction  in  the 
cross  section  v/hich  is  consistent  with  the  Doppler  behavior.  The 
general  trend  in  the  variation  with  time  is  such  that  the  experimental 
cross  sections  hold  up  closer  to  the  peak  for  longer  time  than  the 
synthesized.  This  could  happen  if  the  expansion  phase  of  the  plasma 
frequency  contour  supporting  the  reflection  lasts  longer  than  that 
associated  with  the  Gaussian  model.  The  defocussing  effects  of  the 
underlying  ionization  which  build  up  in  time  are  expected  to  offset 
to  some  extent  the  increase  in  the  cross  section  associated  with  the 
extended  phase  of  expansion  and  thus  can  lead  to  the  observed  slow 
variation.  It  emerges  from  the  comparison  that  the  true  distribution 
of  the  cloud  tends  to  become  flatter  in  the  sense  that  the  density 
falls  off  much  more  gradually  from  the  center  than  that  in  the  adopted 
Gaussian  model  as  the  time  progress  towards  the  signal  drop-out. 
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SECTION  7.  SUMMARY  AND  CONCLUSIONS 


A  three  dimensional  ray  tracing  program  has  been  used  to  syn¬ 
thesize  hf  radar  signatures  from  Gaussian  plasma  clouds.  The  prime 
objective  of  the  investigation  is  to  test  the  validity  of  the  assumed 
Gaussian  model  by  comparing  the  synthesized  signatures  with  the  ex¬ 
perimental  observations.  The  ray  tracing  results  have  been  compared 
also  to  that  computed  for  the  case  of  a  hard  expanding  ellipsoid 
model  to  study  the  effects  of  the  ionization  lying  below  the  level 
of  reflection. 

The  3D  ray  tracing  program  used  in  this  study  has  been  developed 
at  ESSA  by  Dr.  R.M.  Jones  and  is  based  on  a  set  of  six  partial  differ¬ 
ential  equations  similar  to  that  of  Haselgrove  (1954) .  The  program 
can  take  the  electron  collisions  and  the  Earth’ s  magnetic  field  into 
account  and  trace  ray  paths  for  both  ordinary  and  extraordinary  modes. 

The  Gaussian  distribution  of  electron  (ion)  density  in  the  plasma 
cloud  used  for  synthesizing  the  radar  signatures  is  based  on  the  ambi- 
polar  diffusion  model  given  by  Holway  (1965) .  The  Gaussian  parameters 
adopted  for  the  synthesis  are  that  appropriate  to  the  test  Mulberry 
of  Pre-SECEDE  series  conducted  at  White  Sands,  New  Mexico  (Minkoff, 
1970)  . 

The  ray  tracing  method  of  synthesizing  a  radar  signature  from  the 
cloud  involves  first  the  determination  of  the  'homing'  ray  when  the 
transmitter-cloud-receiver  configuration  is  given.  The  Doppler  shift 
of  interest  is  the  one  associated  with  this  ray  and  it  is  recorded. 

The  next  step  is  to  evaluate  the  scattering  cross  section  of  the 
cloud  which  has  oeen  done  by  following  the  technique  described  by 
Croft  (1967)  and  Georges  and  Stephenson  (1969) „  It  involves  trans¬ 
mission  of  a  narrow  bundle  of  rays  and  the  determination  of  the  area 
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of  their  impact  on  the  ground.  The  radar  signature  synthesis  by  the 
hard  ellipsoid  method  involves  the  determination  of  the  specular 
point  from  the  radar.  A  method  to  find  the  specular  point  on  the 
cloud  and  to  evaluate  the  Doppler  shift  and  the  scattering  cross 
section  has  been  described. 

The  experimental  data  for  which  the  ray  tracing  simulation  has 
been  made  was  collected  by  operating  hf  phase  path  sounders  at  10.31 
and  16.078  MHz  during  test  Mulberry  of  Pre-SECEDE  series.  The  data 
has  been  reduced  to  the  form  of  power  spectra  to  measure  the  Doppler 
shifts  and  to  estimate  the  radar  scattering  cross  sections  of  the 
cloud. 

The  following  conclusions  can  be  drawn  from  the  foregoing  study 
of  the  synthesized  and  the  observed  radar  signatures  from  the  plasma 
clouds : 

The  radar  signatures  synthesized  by  ray  tracing  for  the  Gaussian 
plasma  clouds  are  found  to  be  significantly  different  from  those  ob¬ 
served  experimentally.  The  synthesized  Doppler  shifts  are  somewhat 
lower  and  reveal  a  negative  tail  in  contrast  to  the  observed  Doppler 
behavior.  The  scattering  cross  sections,  however,  are  10  to  20  dB 
greater  than  those  derived  from  the  observations.  These  differences 
tend  to  suggest  that  the  electron  density  variation  outward  from  the 
center  of  the  cloud  is  more  gradual  than  that  in  the  adopted  Gaussian 
model. 

The  comparison  between  the  signatures  synthesized  by  the  ray 
tracing  and  the  hard  ellipsoid  model  shows  that  the  effect  of  the 
ionization  below  the  reflection  level  is  to  introduce  a  significant 
positive  shift  in  Doppler  and  a  considerable  reduction  in  the  scattering 
cross  section  by  way  of  defocusing.  The  factor  by  which  the  scattering 
cross  section  is  reduced  is  found  to  vary  from  about  3  at  the  time 


APPENDIX  A.  GAUSS  X  ROUTINE  AND  SAMPLE  RAY  TRACE 

The  subroutine  Gauss  X  listed  here  calculated  the  Appleton  Param¬ 
eter  X  (=  f2/f2)  and  its  time  and  spatial  derivatives  in  geomagnetic 
n 

dipole  coordinate  system.  The  table  that  follows  the  routine  shows 
the  sample  printout  of  the  ray  tracing  program.  The  results  are  for 
the  homing  extraordinary  ray  at  the  operating  frequency  of  10.31  MHz 
and  at  the  time  180  seconds  after  the  release  of  the  cloud.  Figure  A-l 
shows  the  corresponding  ray  plot  in  the  vertical  plane  defined  by  the 
line  joining  the  transmitter  and  the  receiver.  The  portion  of  the 
ray  path  within  the  cloud  where  there  is  significant  bending  is  shown 
alongside  on  an  enlarged  scale. 


1  _  S  u  8  R  Q  U  T.l  N£  -  0  A  U  S3  X _ . 

2  DATA  PI/3. 1M159265/ 

3  REAL  NQiK,  j  l5L*J{LJ3tNUL  _  __ 

*  COMMON  KU)/Wft/!0(  1H)  iqUH»«H00)/XX/XiPXPR(3)  iPXPT 

5  -E.Q-Y  I  VAl_£NcE  <  F  <  3  j  I  PUN.lAJ  L3  Ll J  ITHN  8  rt  (  l  5  )  J  I  (  Re.,  H (  1  9  )  )  i  ( NO  i  * ( 2U  1  ) 

6  1  )  ,  (  r  *  a  (  202  >  >  •(tAu»W<203>>  ttHOtWfZOHn  t  (OR,  *  (  205)  )  ,  (  DZ  » W  l  21)6  )  )  , 

7  2(  K.«H  20  7)1  t  (THn.’M  208  )  3  .  (  PH.O  .  JflLt  2Q-9  LI..MU210)  )  .  (H.  trill  1)J _ 

a  ENTRY  ELECTX 

9  I  F  (  W  (  30  ll,  6  T*  Q«.i  T.Q_  l- _  ... 

|U  W  (  3  0  1 ) » 1 • 

11  kC0S0»hO**2*R • •  0R*T 

12  ZCO»SQRT(hO*«2  +  ‘(.*OZ«T) 

li  CTHH0»S1NiTH0)»S|N(ThN1*CQS(TH0)*.CQ5(THN)»cQSCPHQ»PHN» _ _  . 

1R  STriM0*sQRT  (  1  .UC1-CTHM0**2  ) 

|  5  SPHM0-CQS.I  THUJ  «S.lNJP.«Q-PHN  1VS  T.HMQ 
16  CPHMQ»S0RT  (  1 .0C1-SPHM0**2  ) 

\)  XHO*  (RE  +  H  1  f  STHMOf-CPHMO  . 

1«  YM0»(RE+H)*STHm0*SPHm0 

1-?— .  --7-HL13-LR  E.tilXtX  XHJdJJ _ 

20  Ch1M0*PI/2,-I*aCoS(CthM0J 

21  CCHIM0»CQS(CHIm0J. 

22  SCHIM0«SIn(Ch!m0) 

23  CPCX-CPHM0*CCHf,M0 
2R  SPCX»SPHM0*CCH1M0 

25  CPSX»CPHM0«SChtM0 _ 


81 


26  SPSX«SPHMO*SCH!MO 

2  7  FT?  LL* -EXP .(.-I Z-TAUJ. L/_LIC OAR  CO _ 

2 a  OFOT-l ,/ (TAU*ZCO*RCOSQ  )-FT* ( l ./TAU*M.*DR/RCQSQ*2.*DZ/ZCO*»2) 

29  CONS_T»X/F.**2 _  _ 

30  CNST 1 “NO/p J  ♦♦  t  ,5 

3_1 _ CNST2*C0NsT»CNST  1 _ 


32 

33 

3  H 

3b 

36 

37 

1  CR2»C05(R{2) ) 

SR2*SlNlR<2 J ) 

CR3«C0S<R(3) ) 

SR3-SJN.1RL3)  _ 

XM**R ( l } »SR2*CR3 

YM*R ( 1 ) *SR2*SR3 

38 

ZM*R( 1 ) ♦ C  K  2 

39 

XC“CPCX*  (  XH“XJjgj.tSPCX«  (  YM-YMO  )-SCHLMP*  tZH^ZMO) 

MO 

YC“-SPHMQ* ( XH-xMO  *  +CPHH0# ( YM-YMO ) 

Ml 

ZC»CPSX*  (XM-XHO)  *SPS_X_*lYMrYMP  (Zfl-ZMO  ) 

M2 

GT»EXP(-(xc»»2+Yc*»2)/RC0SQ) 

M3 

ht?exp ( - ( zc/zca ) **2 ) 

MM 

N“CnST 1 *FT*GT*hT 

Mb 

x°const*n 

M  6 

D  (  10)«-N»2**XC/RCOSW 

M  7 

dUH*  -  NA2.1  *jr£/R  cP5_<i_ 

M  8 

D(  )2)«-N*2»*ZC/ZcQ**2 

M  9 

0  1  13)»CPCx*DC in)-SPHM0*0( 11 >*CPSX*D(  12) 

50 

0  (  1M)«SPCX*0<  lnUCPHM0*0(  11  )*SPSX*D(  12) 

bl 

Q  (  1  5  )  a  -  S  C  H  J  HQ  *0  (  LPi+_C  CHiMO  *0(12) 

52 

0  (  l  )»SR2*cR3 

53 

D  (  2  )  *R  (  1  )  * C R 2 •  CR 3. 

5  M 

D  (  3  )  *-R ( 1 )*SR2*SR3 

5b 

0(M)-SR2*SR3  . 

56 

0  (  5  )  "R ( 1 )*CR2*SR3 

57 

0  (  6  )  *R  (  1  )  •  S R 2 •  c R 3 

S8 

0(7)-CR2 

59 

D  (  8 ) *-R ( 1 ) *SR2 

D  (  9  )  *0  « 

61 

P  X  P  R  (  l  )  *C*ON  S  T  •  (  0  (  1  3  )*  0  (  l  )  *  D  (  l  H  )  *  0  (  M  )  ♦  D  (  15)tD(7)  ) 

62 

PXPR  (  2  L*?_OJMST*_(Ji  ( J  3J_*  p_(  2)±  CM  1  M  )  •  0  (5  ) +0  (  l  5  )  *0  (  8  )  ) 

63 

PXPR(3)»C0NST*(0( 13) •0C3)*0* H)*0(6)  ) 

6  M 

0  (  2  )  "GT*M  «  *QR* ( (XC**2*YC**2)/RC0SQ*»2) 

65  ' 

D(3)»HT*M,*DZ*(ZC/ZC0**2)**2 

66 

p xpt®cnst 2 •  Ui?ht*pfj>i*ft.1HT *P±2! *fT*gt*o ( 3)  > 

67 

RETURN 

68 

END 

82 


1 


to 


Ic  c 

IW  ICO 

_J  v'o  o 

ia  st©  o 

ft°l  *! 


sec  e!c  c 
O  Op  OlO  c 
O  O  O  O'O  c 

s  o;o  0,0  o 


I 


» 


1  I 

!i  ei 


to.  cd!  • 

|gci 

lm  i 
i«  I 


l 


ia  !c  o 

lm  j  £ 
lx  t 
la  I 

IX  io  £ 

*-  ,o  m 


I 
i 

i  i 

kc  £*£  S C*£  £l 
b  mm  mm  to 

Ln  .o 


to  .w  — 
|a  too 

!i5i"4 

!5  ! 
i  2  *0  m 
oio  o 

,UJ  I-1— 

(-1  <11 

££b 

iS5i°. 

'  t!1 

'  Si 

'  * 
Zl«*  Of  £ 

2>c  c! 

5J  t  - 
>1  r- r- 
jto  oM 

uiss«« 

Ik  |f»  ft 


I  »' 

I  ' 

Jc  £  :£  vOlvO  £ 

k>  m  m  mm  m 
l£  £>£  £l£ 

^  ft,ft 


»  I 


lo  oio  o  o  o 


I  HI  1(1 

bob  ob o 
poo  op  o 


.  .:lo  u|  • 
F2lf>°!  * 

S'  ' 

r—  *-*l  |  O 

p  ujla  ol  o 

*iza  : 

H  I 
I  I 
J  io  ft 

lid  IC  4 

B 

•  i  * 
i 

I  'o  3 

I  B3 

•UION 

P‘!*o 

U  {  0 
I  I  "* 
I  I  , 
I  3 


I 


I  I 


I  to  o 


•  I .  .1  • 
I 


,£  Wm  mb  m 

m  >o  to  ft, is  to 

£  £'—  —  o  £ 
!r»  ftto  ©]g  a 

(*b  gift  g'ft  a 


—  m,to  ft.o  ft 
r-  ft  r»  ft'a  x 
££fo  £|£  p 

to  £  in  n'n  £ 

-  ftjft  r-jft 

^  ft  r»  ftg-  ft 
oo.g  a, a  c, 

ft  ft 'ft  ft’ ft  ft 

U  Xl 


pop  0(0  o 


[o  OIO  0,0  o 
p  o  .o  0.0  o 
o  o!o  oio  o 


U  »  }0  g1*  _ 
cr  -|£  fttm  - 

p  ft  to  win  ft| 

-  —  ioj  ftift  ft 

t  1\T  t|t  3 


_  £'n  n 

to ft  |c  £3  m 
pf'iJ  -  o  £ 
OK)  lO'tV  <V 


Irj  >0  j  f* 
0  0,0  0(0  9 

“r-r- 


p  oio  o  o  o 

r  •  I  II  I 


c  c'A  m#  r* 
[c  o,c  c.vr  o. 
[O  o'o  o'o  ~\ 
o  0,0  o'o  o| 
*1  *  *1 
i  > 


c  o  aV  r» 
p  -:C  g.r*  ^ 

m  tr'j  -  ©  - 


,  _  JV  <  fV 

fc?<CSCi? 

U  -  NM« 


I 


<m  «>•  f**a 
P  *o:o 
<c  r-,—  r*  ^  c 
e  *o‘c  a'j  c 

a  0.^4  m[o  <o 
c  a!c  cla  a 


•£>  i/>  «c  -«  ■'r^r*r^aa 
O  c'3  I/I.X  o  £>  C  >C  'CSC  X, 
o  vT’.c  » »r  tfi  n»r*n  kr^r  irj 

»  »  j  o  i/>  /i  ft  in  ^  in 
^  o  p  g  p  c  c  v  o  a  ,o  o 


p  o.b  x>!JO  >0  n  -ftp  £ 

|j  3*2  22  22  22  2 

*  r-;o  o  o  oio  o  o  D  o  o 


a  <P«  -,tr  ec 
x  «lo  a 
■o  «ojx  r*  n  « 


,  f-jo  9', 

o-  &‘o  o‘a  a, 
icv  <V.f\J  MlftJ  <v| 

n  Ji'n  .ft!«n  n 

r*  -  ' 

h  «;r*  a. 

O  OP  0J,«O  H 

o  o'o  o!<vxJ 

I*  I*  I  .  1 1 


O  op  Op  X, 
o  o'o  o  o  o 
o  so  o  o  o 


o  Ol-* 
Mli/ 
p  Pr- 
P  <Ml^* 

|*C  «ft(© 

t  a^a 

I 

^  *>l» 


-4I  r 
r;o 
/»  o  _. 

OJifV  x] 

tin  n 
*1* 

I 

J1!  J-  n| 

rir>.  © 
3>  o| 


O  Op  f 


to  an bp « 
K  C'o  oC  . 


p  «np  oxx 


I  lit  I  >  I  I 


O  O  O  Or©  O 


I  I  fl  I  <1  I 

b  o  b  o  o  © 

p  op  op  - 


p  a>r* 
N-OOMJ, 

|o  xi.t  o  >0  n 

X  2'A  <c!«  <C 
> 

p  xi-*  a 
o  pm  in,r*  x 
0  *o  »o  m  m  m1 


k  -ix  x‘m  x 

fN.  — 

j  t'o  x  x  t 

1  np  m!x  g 
1 ‘i  - 
,  *H  » 
lo  oj—  op  ml 
p—  C 
3  o'—  fti'm  m 


g  n  op  ft. 
,m  mp  x'm  x 
t  t  t  -v  n  t1 
I-  aim  x  ^  x 

C  *  s'  ■pT  s' 
*  t‘t  »*t  t 


c  o'm  t  vn 

T  y|,M  O'?  C 

g  -.0  »p  m 
g  xt<v^--t  c 


I  I  I  I  >1  I 


£  X  £  t*  tr  >r 

3  t '»  £  £  £ 
c  £  ;£  —  k3  9 
?  3  3  «T  £  £ 


3  OOOOO 
O  O'C  c  c  c, 
ooioopo; 


n  o  o  o  ,n  ift 


O  £  -O  £ 
p  mm mm m 
/-*—  ■-*—  *£  £ 
in  n  n  o  tv  rw 


P  m  m  3  *n  n 
|<v  <v  tv  rv  i<v 


p  o  sV  vo  x> , 
Jo  —  n  tv  nn' 
i3  3 13  nm  n 

In  nin  n^  n 
►-•  —P  r»  —  — 

kc  tv<v  a  o  o 


--  —  a\r  n 
3  3  3  m  o  c 
£  £,£  <V  •£  «£■ 

««,«irnm| 
—  —  n  r* 

•  11  I  ’I  II 
-*  —  on  n 
£  £i£  £  V  U* 
In  /ip  in n  n 


bmmnoo 

bft.xv£- 


•  •••  *  • 

K  ft.  3  n  n 
to  fti'o  g  o  £: 
3  n'£  3  ■ 


,t  oj>j  1  -* 
x  o  m  £  n  n 
loc«  \  t  t 
T  O  —  M  ft,  \ 

*  x  r  tc  j? 

»•»  if»  " 

1 

r*-  t  m  o  o 
9  CP  Op  V 
o  —  O  —  -*’ 
tc  n>- 


r  ©U  vm  m 
[b  g  a  g  g  g 

r 


pggg-*- 
In  m  m  m  c  c 
C  f  0  £*•  r. 
p  <vrv  <vrg  tv 


poo  V»3  t 

to  cx  \i o  o 
b  00  o  £  9 

p  ft»  ■%  ft»«-*  — 


m mm mm  m 
b-  gg  gg  g 

-W,  -.1 


O  > 

*  i 


I  S  <  t<  t 


I  I  I  I  I  I 


>  £  m  r  c 

x  r-  r*  c  3  o 

3  r-  £  fti  3 

£  c  o  n  c  g 


ft  mk  m  3  » 
o  —m  o  £  ft 
£  tn  9  mm  ft 

9  g  g  g  g  g 


«  ft  o 
a  f>  x 
SCO 

g  g  g 


poo  op  o, 
c  c’c  o’c  o 
O  0,0  OIO  c 

o  ob  000 


«  m  r*  ft  r  c 
to  £  — •  ft*  a  3 
fo  3  m  <  3  «r 
C  3  m  ft-  rg 

bfvm  m*3  £ 


£  £  £  £,£  £ 
m  m  m  m  m  ml 
£  —j—  —  —  — 
NbbOOO 


£  £  «£  £  £  £ 

n 22  222 

jo  OO  tftOO 


v>ft»  a  i/vg  m, 
—  —  —  01,01  m 
ft  ft  Ol  ft  ft  ft 


tv  bib  3  m  ft 


KV  ftlft  ft  ft  ft 

b  o  o  o  o  o 


n  *  «  nm 
con  mm  £ 
c-o  m  —  3 


ban  —  o  r- 
ft  £  a  c  m  m 

to  — o  00  g 


m  x  g  o 
g*  ftf*  «.«  « 1 
11  11  1  1  1  1 
nr-m  n  n  g, 
in  m  —  o  g  £ 
n  bib  3  o»  — 


b  3  —  3»C ft 

lo  or-  ir  m  r* 
Jo  og  g  a  x 


omr«  m  mn 
o  c-  g  o  3, 
—  ft£  —  x  «n 


tv  —  m  n  n  ft 
'  £  g  g  o-  x 
-o  g  n- ft 


g^  coo 

Xf-  ft  3 
x  a  m  r  3 
ft  ft  m  m'ft  o 


3  om  0  3  0 
st—  cm  — 
c  c  —  *»  n  9 
-  gr-  ft  3  x 


\  t  ft  m  -g  n 


C  CO 

000 
pop 


kv  ft  ft  ft  ft  ft 


ft  ift'  o 
c  —  — ' 
9  g  o, 

O  Of  oj 
3  3 

ft  ft  cvl 


£  £  £ 
1  mm 

-*  — 1  — 
000 


£cc g  p  ft 
to  mm  m  3  3 , 
kv  o.ft  ft 01  ft 


to  O;  3 

L5  n  n, 

K.ft<ft 


b  «  ft  x  •£  x 
tft  g  3  —  o  o 
U>oo  3  O  £ 


k>  g  3  o£  3 
‘c  m  -  n  c  ft 
C*-  k  r.  g 


0  Oft  >  o  ft 

it*,  r-  £  o,  g  m 
K  m  -  n  r  — 
r*  £3  mrg  o 


c  g»c  3  v ft 
»  ft  in  ft*-*  n 

t  om  CTO 
3  Ch  ftft*  X 


s , 


In  n  m  t  t  r 


€  *'<  XI 

£  px  £ 


ft  ftft  ft 
I  l.l  1 

I  I 


040  a 
o  0,0  o 
•5  •  •( 


B 


._  o  £ 

bS5| 


Fjr-I 


.t  x  ft 
Jn  g  - 
t  g  n 


1  l 

a  g’an  £i 

o  OJO  c* 

g  0  £  £ 
o  o'm  n 

«,  •  •> 

hO 

ft  ftIC  <3 
ft  orm  m 


£  £;ft  Ol 

I2SS2 

O  Ojft  ft 

m  mm  f 


£  toft  ft 
ft  T 

I  I 

ft  0/01  oi 
o  00  a 


o  00  c 


o  ob  ot 
o  op  o 


o  oio  o 
®  fl®  «? 
g  g*g  g 
«  «r  g  c. 
o  Po  o 

-  •  r 

m  m  £  £ 
n  ton  to 
£  tox  <r 


C  «,£  kO 
r.  p.  •» 
mm  -  -I 


g  t  o  tt 

—  *nm  m 
c  t  »  t 
g  g'oj  ol 


m  m  o  o 
£  £;©  o 
0  o  o  p 


NOT  REPRODUCIBLE 


83 


lv  t  Real  Part  (n  )j"  lj  Whare  V  la  th«  Magnitude  of  th#  wave  normal  vector  and  n  ia  the  coaplax  phase 
refractive  index.  Thia  quantity  would  be  zero  if  there  were  no  errors  in  the  numerical  integration. 


Figure  A-l  A  Sample  Ray  Trace  Projection  Plot  In  A  Vertical  Plane 
Defined  By  The  Line  Joining  The  Transmitter  And  The  Re¬ 
ceiver.  The  Boxed  Segment  Where  The  Ray  Suffers  Signi¬ 
ficant  Bending  Is  Shown  On  An  Enlarged  Scale 


APPENDIX  B.  A  VECTOR  METHOD  TO  CALCULATE  CLOUD  TRAJECTORY 

The  vector  method  described  below  determines  the  coordinates  of 
the  cloud  using  the  line  of  sight  data  given  in  azimuth  and  elevation 
angles  from  a  pair  of  observing  stations.  The  cloud  position  is  de¬ 
fined  ideally  by  the  intersection  point  of  the  lines  of  sight  from 
the  two  sites.  In  practice,  however,  slight  errors  in  sighting  cause 
the  two  lines  to  miss  each  other  and  in  this  case  the  midpoint  of  the 
common  perpendicular  to  the  two  lines  of  sight  is  taken  to  be  the  best 
estimate  for  the  cloud  position. 

Let  I  and  J  be  the  unit  vectors  along  the  two  lines  of  sight  and 
N  be  the  unit  vector  along  the  common  perpendicular  to  the  two  lines 
such  that 


N 


I  x  J 
I  x  J 


(B-l) 


and  R2  represent  the  vectors  to  the  two  observing  sites  from  the 
center  of  the  earth;  X  denote  the  vector  to  the  cloud  center.  and 
B2  are  the  scalar  distances  from  the  two  sites  to  the  points  joined 
by  the  common  perpendicular.  D  is  the  distance  between  the  two  lines 
along  the  perpendicular.  Then  from  the  geometry  shown  in  Figure  B-l 
we  have; 


^  +  B  7  +  j  N  **  X  (B-2 ) 

*2  +  B2  ^  f  ^  =  *  (B-3) 

Taking  the  scalar  product  of  Equations  B-2  and  B-3  with  N  and  sub¬ 
tracting  one  from  the  other  would  result  in  the  relation 

(*1  "  V  * (  N  )  =  D  (B-4) 
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Figure  B-l  Geometry  for  the  Vector  Determination 
of  the  Cloud  Coordinates 

It  is  required  now  to  evaluate  either  or  B 2  to  obtain  the  posi- 
tion  vector  X  of  the  cloud.  This  can  be  done,  for  example,  by 
vector  multiplication  of  Equations  B-2  and  B-3  with  the  unit  vector 
I  which  yields  the  following  expression  for  B^  on  making  use  of 
Equation  B-l. 


r  x  J  N  =  Tx  [IT  -  R^]  -  D  [n  x  i] 


(B-5) 


Using  Equations  B-4  and  B-5,  the  vector  X  can  be  expressed  in  terms 
of  I,  J,  N,  and  which  can  be  evaluated  from  the  known  quanti¬ 
ties,  the  coordinates  of  the  observing  sites  and  the  look  angles 
to  the  cloud,  as  follows: 

Let  (•')^,  <f>  )  and  ke  the  cO“latit.ude  and  east  longitude 

of  the  two  stations  in  the  geographic  co-ordinate  system.  The 
azimuth  and  elevation  angles  to  the  cloud  from  the  two  stations  are 
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t 


denoted  as  (a  ,  ft  )  and  (a0»  P  ),  respectively.  Define  two  Cartesian 
coordinate  systems,  ( x,y,z )  an  Earth  centered  system  and  (x^ ,  y  ,  2^) 
with  its  origin  at  one  of  the  observing  stations  (say  site  1)  as 
shown  in  Figure  B-2.  (x,  y,  2)  is  such  that  2  is  pointed  toward  north 

and  (x,  2)  defines  the  reference  meridian  plane.  (x^,  y1 ,  2^)  is 
tilted  with  respect  to  (x,  y,  2)  by  0X  in  longitude  and  (90  -  0^  in 
latitude.  (i,  j,  £)  and  (1^,  3  ,  ic  )  are  the  unit  vectors  along 
(x,  y,  2)  and  (x^,  y  ,  2^),  respectively.  They  are  related  by  a  trans¬ 
formation  matrix  as: 


H-> 

H 

_ 1 

sin 

COS  0 

sin  0^ 

sin  0 

cos 

ei" 

A 

31 

= 

-sin 

*1 

cos  0^ 

0 

-cos 

01 

COS  0^ 

-cos  0^ 

sin  0^ 

sin 

Z 


1 


A 

D 


ic 


(B-6) 


Figure  B-2  The  coordinate  systems  involved  in  the  transformation 

of  various  vectors . (x,  y,  2)  is  an  Earth  centered  system 
and  (x^,  y^,  z^)  is  centered  at  an  observation  site. 
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The  unit  vector  I  along  the  line  of  sight  from  station  1  is  given  in 
(x^,  ^l'  Z1 ^  s^stem  terms  t^*,e  1°°^  angles  as: 

I  =  sin  J*  i  +  cos  sin  3^  +  cos  0^  cos  (B— 7 ) 

substituting  Equation  B-6  into  B-7  leads  to: 


=  (sin  3^  sin  0^  cos  0^  -  cos  |3  sin  sin  0^ 


-  cos  3  cos  «  cos  0^  cos  0^  )  1 
+(sin  3^  sin  0^  sin  0^  +  cos  3^  sin  cos  0^ 


(B-8) 


-  cos  3^  cos  «  cos  0^  sin  0^)  j 

i  A 

+  (sin  m ^  cos  0^  +  cos  3^  cos  a  sin  0^  )  k: 


}  ^  A 

=  a^i  +  b1  ]  +  c^k 


J,  the  unit  vector  from  the  second  site,  can  be  similarly  expressed; 
replacing  the  subscript  1  in  Equation  B-8  by  2  one  gets: 


J  =  a  i  +  b~  j  +  eft 


(B-9) 


From  Equation  B-l 


A  ■!. 

i  3 


N  =  j  a1  bx  Cj 


a  b  c 
2  2  2 


I  x  J 


Nx  i  +  N2  D  +  N  k 


(B-10) 


The  vectors  R^  and  R^  are  given  in  terms  of  the  coordinates  of 
the  stations  as 


1=  R  sin  cos  0.  l  +  sin  F>.  sin  0,3  +  cos  a.  k  I  (B-ll) 

1  E  [  1  1  1  1  1  j 
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R2  =  Rg  (sin  cos  02  i  +  sin  @2  sin  <p2  j  +  cos  ©2  ft  )  (B-12) 

using  Equations  B-4,  B-5,  and  B-8  to  B-12,  the  position  vector  of 
the  cloud  can  be  expressed  in  the  Earth  centered  coordinate  system 
in  the  form: 

X  =  U  1  +  V  ]  +  W  (c  (B-13) 

Where  U,  V  and  W  are  functions  of  the  parameters  9  ,  9  ,  0  ,  0«,  , 

J.  6  X  ^ 

a 2<  ^ i  an<3  an<^  can  caicuiated.  If  h,  9,  and  0  denote  the  al¬ 
titude,  latitude  and  longitude  of  the  cloud  position,  they  are  given 
as : 

h  =  (U2  +  V2  +  W2  )  1/2  -  Rg 

0  =  si,."1  (W(02  +  v2  +  wV1/2) 

0  =  tan-1  (v/u) 

The  method  described  here  has  been  used  to  calculate  the  trajectory 
of  the  cloud  shown  in  Figure  3-2. 
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